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SECTION  I 
INTRODUCTION 


This  report  provides  user  instructions  for  the  1987  version  of  the  three-dimensional 
EPIC-3  code. 

The  first  documented  version  of  the  EPIC-3  code  was  developed  for  the  Ballistic 
Research  Laboratory  (BRL)  in  1977  (Reference  1),  the  second  version  was  developed  for 
the  Air  Force  Armament  Laboratory  in  1978  (Reference  2),  and  the  third  version  was 
developed  for  the  BRL  in  1980. 

The  1987  version,  described  in  this  report,  incorporates  a  NABOR  node  option  for  vari¬ 
able  nodal  connectivity.  Other  new  features  are  improved  material  models  and  a  data 
library,  which  are  consistent  with  the  1986  version  of  the  EPIC-2  code  (Reference  3); 
eroding  interfaces  for  projectile  penetration  into  thick  targets;  and  a  target  add  option. 
Numerous  other  changes  have  also  been  incorporated  since  1980. 

One  of  the  distinguishing  features  of  the  EPIC-3  code  is  that  it  is  primarily  based  on  a 
tetrahedral  element  formulation.  This  allows  problems  with  severe  distortions  to  be 
run  without  rezoning.  This  is  a  very  desirable  feature  for  many  intense  impulsive  load¬ 
ing  problems  involving  high  velocity  impact  or  explosive  detonation. 

The  1987  version  of  the  EPIC-3  code,  as  well  as  the  user  instructions  in  this  report,  have 
been  made  as  similar  as  possible  to  the  corresponding  1986  version  of  the  EPIC-2  code 
(Reference  3). 


SECTION  II 

USER  INSTRUCTIONS 


This  section  provides  user  instructions  for  the  EPIC-3  code.  The  code  consists  of  a 
Preprocessor,  Main  Routine,  and  Postprocessor  for  state  and  time  plots.  The  formula¬ 
tion  is  not  provided  here;  however,  most  of  the  equations  are  identical  to  those  of 
Reference  2.  Some  comments  concerning  the  formulation  of  the  current  version  of  the 
code  follow: 

•  The  material  strength  models  for  solids  and  crushable  materials  are  given  in 
References  4  and  5. 

•  The  fracture  model  for  solids  is  given  in  Reference  6. 

•  The  average  pressure  formulation  is  given  in  References  7  and  8. 

•  The  sliding  interface  formulation  is  given  in  References  9  and  10. 

•  The  eroding  interface  formulation  is  given  in  Reference  8. 

•  Examples  of  anisotropic  computations,  performed  with  the  1980  version  of  the 
EPIC-3  code,  are  given  in  Reference  7. 

•  The  three-dimensional  NABOR  node  formulation  will  appear  in  the  final  report 
for  this  contract.  The  NABOR  formulation  for  plane  strain  is  given  in  Reference 
11. 

•  A  characteristic  of  tetrahedral  elements  is  that  they  sometimes  produce  a  pres¬ 
sure  field  that  oscillates  spatially.  This  can  cause  inaccuracies  for  cases  where  the 
material  strength  or  fracture  is  dependent  on  the  pressure  in  the  element.  This  is 
overcome  by  computing  a  nodal  pressure  which  is  the  average  of  the  element  pres¬ 
sures  for  all  elements  that  contain  a  specified  node.  The  nodal  pressures  do  not 
oscillate  spatially  but  rather  form  a  smooth  pressure  field.  A  smoothed  element 
pressure  is  then  defined  as  the  average  of  the  four  associated  nodal  pressures.  This 
smoothed  pressure  is  used  for  the  strength  and  fracture  models,  which  are  pres¬ 
sure  dependent. 

A  description  of  input  data  for  the  EPIC-3  code  is  given  in  Figures  1  through  20. 
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isssssmsssss^^ 


SPHERE  ELEMENT  OESCRFTION  CARO  (26.5X.26) 


SPHERE  ELEMENT  MATERIAL  CARO  FOR  MATL.  0  (166) 


25321 


NABOR  FILLED  PLATE  ELEMENT  IDENTIFICATION  CARD  (6) 


issssssss^^ 


NABOR  FILLED  PLATE  ELEMENT  DESCRIPTION  CARO  (76) 
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Figure  4.  Element  Input  Data 


DESCRIPTION  CARD  (A80) 

I  DESCRIPTION  OF  PROBLEM 


l _ _ _ _ 

IDENTIFICATION  CARO  (SIS,  5X.  5FI0.0) 

l  gall  320  ■ 1 ; '  :’ci  ■ 

VFRACT 

|  PMAX  |  VNREF  |  ARMAX  |  ASERR  | 

NTEGRATION  TME  NC  RE  ME  NT  CARD  (6F10.0) 

|  DTMAX  1  OTMM  I  SSF  I 

TMAX 

PRESSURE  CARDS  FOR  IPRES  -  2  -  AS  REQUIRED  (SIS,  FI 0.0) 

PRES 

|  BLANK  CARD  ^  NCLUOED  ONLY  FOR  IPRES  -  2 
TME  -  PRESSURE  CAROS  FOR  IPRES  -  2  -  AS  REQUIRED  (2F10.0) 


PTME 


P<T) 


I  BLANK  CARD  ^  INCLUDED  ONLY  FOR  IPRES  -  2 
CHANGE  DROP  CARO  FOR  CHANG -1QR  2  (F10.0,10B) 

|  TDROP  f/OOE  fopoEf  ELE  |  PELE  INSLID  H«-OT>.PLOT  j  NFAIL  |  NAB  |  PNAB 
DESIGNATED  ELEMENT  FAILURE  CARO  -  FOR  NFAIL  >  0  (1615) 

I  EF1  I  EFZ  I  I  I  I  I  I  I  1  EFN  i  I  I  ii  1  - J 

CHANGE  ADO  CARO  FOR  CHANG -2  OR  3  (F10.0,  It  g) 

|  TAPP  |NSLIO  P(YZRIG|  SPLIT  |  N1R  |  NNR  |  NIL  |  NOR  |  NR  [  NPL  |  NRSTl PRINT 

TARGET  SCALE /SHIFT /ROTATE  CARD  FOR  CHANG.  2  OR  3  (7F10.0)  _ _ 

|  XSCALE  |  Y SCALE  |  ZSCALE  |  XSHIFT  ]  ZSHIFT  |  ROTATE  |  SLANT  bvSSSSSSSS^l 

ADDITIONAL  TARGET  NODE  DATA  FOR  CHANG- 2  OR  3 -AS  REQUIRED  ' 

[  BLANK  CARD  ^  ENOS  NODE  DATA  FOR  CHANG  -  2  OR  3 

ADDITIONAL  TARGET  ELEMENT  DATA  FOR  CHANG  -  2  OR  3  -  AS  REQUIRED 
I  BLANK  CARD  ^  ENDS  ELEMENT  DATA  FOR  CHANG  -  2  OR  3 

NSLtO  DA™  -  AS  REQUIRED  FOR  CHANG  .  2  OR  3\  g^r-  FORMAT  AS  USED  IN  PREPROCESSOR 
NRST  DATA -AS  REQUIRED  FOR  CHANG- 2  OR  3  J 

TARGET  VELOCITY  CARD  FOR  CHANG  -  2  OR  3  (30 X,  3F10.0) 

■***  i  i 

|  BLANK  CARD  ^  ENDS  DATA  FOR  CHANG  -  2  OR  3 


Figure  5a.  Main  Routine  input  Data 
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PLOT  CARO  (515, 5X,  4F10.0) 

[CYCLE  |  SYS  |  LOAD  ftPLOT  |  LPLOT]^sS3  OT  SYS  |  TSYS  |  PTnOoE  |  TNOOE  * 


DESIGNATED  NODES  CARDS  -  AS  REQUIRED  (166) 


IN1  lN2  I  ~T 


1  11  1 


DESIGNATED  ELEMENTS  CAROS  -  AS  REQUIRED  (166) 


E1  I  E2  I  |  |  |  |  |  T 


EN  |  1  |  1 


DATA  OUTPUT  CARDS  -  AS  REQUIRED  (2P10.0, 36) 


TME  |  ECHECK  [LOAD  |SAV£  (YPRKTE^^^^^ 


BLANK  CARO 


Figure  5b.  Main  Routine  Input  Data 
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GEOMETRY  PLOT  CABO  (215,  FIO.O. 212, 4X,  212, 18X,  A30) 


TIME 
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PLOT  UMfTS  CABO  FOB  AXES  -  2  (6F10.0, 215) 


XMW 

YMAX 

YMW  1 
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ZMIN 
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vELoerry  vector  plot  cabd  <215,  F10.0. 212. 2X.  212.  iox,  fio.o.  aso) 
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PLOT  UMfTS  CABO  FOB  AXES  -  2  (6F10.0) 
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30  PEBSPECTIVE  CABO  FOR  VIEW  -  4  (6F10.0.) 
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CONTOUR  PLOT  CABO  (215,  FIO.O,  212,  2X,  412,  I6X,  A30) 


YMAX 

YMW 
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ZMIN 

CONTOUR  SPECIFICATION  CARO  (8F10.0) 


VAR  (1)  |  VAR  (2) 


VAR(NLWE) 


PRESSURE 

VON  MISES  STRESS 

EQUIVALENT  PLASTIC  STRAW 

DAMAGE 

TB4PERATVJRE 

PLASTIC  WORK  PER  NfTlAL  VOLUME 
WTERNAL  ENERGY  PER  WITIAL  VOLUME 
LOG  (10)  STRAW  RATE 
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PLOT  UMfTS  CABO  FOR  AXES  -  2  (7F10  0) 
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Figure  6.  Postprocessor  Input  Data  for  State  Plots 
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SYSTEM  PLOT  CAROS  -  A3  REQLHRED  (215, 5X,  F5A  4F10.Q,  A2D) 

|  TYP6|  AXESP^s^SCALEI  TMAX  |  TWIN  [  VMAX  |  VMW  | 


PTIERNAL  LOADS  PLOT  CAROS  ■  AS  REQLHRED  (3B.  F5.0, 4F10.0.  A20) 


TMAX 

~|  TMN  (  VMAX  |  VMW  | 

- mj 

WOMDUAL  NODE  PLOT  CARDS  ■ 

AS  REQUIRED  (316,  F5.0,  4f  10.0,  A20) 

TMAX 

|  TMW  |  VMAX  |  VMW  | 

TITLE 

NCXVBUAL  ELEMENT  PLOT  CAROS  -  AS  REOUWED  (3»,  F5.0, 4F10.0,  A20) 
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•ETERNAL  LOADS:  AXIAL  (27).  SHEAR  (28),  BENDMG  (29) 

XYZ  NOOE  POSITIONS 

XYZ  NOOE  VELOCmES 

XYZ  NOOE  ACCELERATIONS 
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DAMAGE 

RATIO:  MEAN  STRESS  /  VON  MISES  STRESS 
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Figure  7.  Postprocessor  Input  Data  for  Time  Plots 
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Figure  8.  Nodal  Spacing  for  Various  Expansion  Factors 
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Figure  10.  Node/Element  Input  Data  Example 
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Figure  1 1 .  Rod  Shape  Geometry 
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Figure  16.  Summary  of  Nodes 
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Figure  18.  Pressure  Model  for  Crushable  Solids.  Specific  Data 
Shown  are  for  Concrete 
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Figure  19.  Euler  Angles  for  Principal  Material  Axes  X',  Y,' Z' 
Relative  to  Initial  System  Axes  X,  Y,  Z 


Figure  20.  Sign  Convention  for  Internal  Loads 


1.  INPUT  DATA  FOR  THE  PREPROCESSOR 


The  functions  of  the  Preprocessor  are  to  define  the  initial  geometry  and  velocity  condi¬ 
tions,  and  to  determine  the  memory  size  requirements  for  the  Main  Routine.  The 
descriptions  which  follow  are  for  the  data  in  Figure  1.  Consistent  units  must  be  used, 
and  the  unit  of  time  must  be  seconds. 

It  is  possible  to  interject  user  comments  into  the  data  by  use  of  a  $  character.  If  the  $  is  in 
the  first  column  of  the  card,  that  entire  card  is  ignored  as  input  data.  If  the  $  is  beyond 
the  first  column  in  the  card,  then  the  $  and  all  data  to  the  right  of  the  $  are  ignored. 

Description  Card  (A80)  —  A  description  of  the  problem  provided  by  the  user. 

Miscellaneous  Card  (1415)  - 


CASE 

PRINT 


SAVE 

NSLID 


XYZRIG 


=  Case  number  for  run  identification. 


0  will  not  print  individual  data  for  each  node  and  element. 

<(  1  will  print  individual  node  and  element  data,  and  complete  slid¬ 
ing  interface  data. 

2  will  print  node  data  for  nodes  with  y  =  0  and  element  data  for 
elements  with  one  face  on  the  y  =  0  plane. 


f'O  will  not  write  results  on  restart  file. 
Ll  will  write  on  restart  file. 


=  Number  of  sliding  surfaces. 

'100  gives  a  rigid  frictionless  surface  on  the  positive  side  of  the 
plane  described  by  x  =*  0.  If  the  equations  ofmotioncauseanodeto 
have  a  negative  x  coordinate,  then  the  x  coordinate  is  set  to  zero 
and  the  x  velocity  is  also  set  to  zero.  NABOR  nodes  with  an  x  coor¬ 
dinate  less  than  one  radius  have  the  x  coordinate  set  to  one  radius 
and  the  x  velocity  set  to  zero. 


10  gives  a  rigid  frictionless  surface  at  y  =  0. 


1  gives  a  rigid  frictionless  surface  at  z  =  0. 

0  or  blank  will  ignore  this  option. 

These  options  can  be  combined  by  addition.  For  example,  111 
Vgives  rigid  frictionless  surfaces  on  all  three  major  planes. 


"0  will  perform  the  sliding  surface  computations  after  the  updated 
velocities  and  displacements  are  determined  from  the  usual  equa¬ 
tions  of  motion.  Contact  is  established  as  long  as  the  slave  node 
interferes  with  the  master  surface  before  the  velocities  and  dis¬ 
placements  are  adjusted  in  the  sliding  surface  routines.  This  op¬ 
tion  is  the  most  reliable  and  should  be  used  for  complicated  slid¬ 
ing  surfaces,  which  include  double  pass  and  intersecting  sliding 
surfaces.  It  must  be  used  for  the  eroding  interface  option.  For  rela¬ 
tively  low  velocity  impact  problems,  where  there  is  limited  defor¬ 
mation,  this  option  can  introduce  significant  errors  in  the  form  of 
excessive  deformation  and  internal  energy.  A  double  pass  sliding 
surface  is  a  pair  of  sliding  surfaces  that  have  the  role  of  Master 
and  Slave  interchanged.  The  two  sliding  surface  passes,  one  from 
each  sliding  surface,  prevent  either  surface  from  penetrating  the 
other. 

=< 

1  will  perform  the  sliding  surface  computations  after  the  updated 
velocities  are  determined  but  before  the  updated  displacements 
are  determined.  This  should  be  used  if  the  sliding  surfaces  do  not 
require  a  double  pass  or  the  eroding  interface  option.  Contact  is 
first  established  when  the  slave  node  interferes  with  the  master 
surface.  Thereafter,  a  slave  node  is  considered  to  be  in  contact  un¬ 
til  the  preadjusted  normal  velocities  between  the  two  surfaces  are 
separating  rather  than  closing.  This  approach  minimizes  the  dis¬ 
tance  the  slave  node  is  moved  to  place  it  on  the  master  surface  and 
is  therefore  more  accurate.  If  there  are  no  sliding  surfaces,  either 
..option  can  be  used  (SPLIT  *  0  or  SPLIT  =  1). 

=  The  first  node  in  a  series  of  consecutive  nodes,  which  act  as  a  rigid 
body.  This  option  automatically  restrains  the  nodes  in  the  y  direc¬ 
tion  and  leaves  them  free  to  move  in  the  x  and  z  directions.  The 
rigid  body  nodes  cannot  be  restrained  in  either  the  x  or  z  direc¬ 
tions  and  they  cannot  be  designated  as  slave  nodes  on  sliding  sur¬ 
faces. 

=  The  last  rigid  body  node.  All  notes  from  N 1R  to  NNR  are  treated 
as  a  rigid  body.  This  node  number  can  control  the  required 
bandwidth  since  the  minimum  number  of  nodes  in  memory  must 
be  at  least  equal  to  NNR 

=  The  node  number  of  the  innermost  node  at  the  top  free  end  of  a 
cylindrical  projectile.  This  is  used  if  it  is  desired  to  compute  inter¬ 
nal  loads  in  a  slender  projectile.  If  this  option  is  used,  the  nodal 
geometry  must  be  consistent  with  the  arrangement  of  nodes 
generated  by  the  rod  node  generator  described  later.  Leave  blank 
if  not  used. 

=  The  outer  ring  number  of  the  rod  if  internal  loads  are  requested. 


NIR 


=  The  inner  ring  number  of  the  rod. 

-  The  number  of  planes  of  nodes  in  the  rod  for  the  internal  loads. 


NPL 

UNIT 

NRST 


'Q  indicates  the  constants  in  the  material  library  have  English 
units  (pound/inch/second/degree  Fahrenheit). 

=< 

1  indicates  the  constants  in  the  material  library  are  converted  to 
^Standard  International  (SI)  units. 

=  Number  of  groups  of  nodes  to  have  restraints  redefined. 


Material  Data  Cards  -  Material  data  can  be  defined  by  the  user  or  taken  from  the 
material  data  library.  Specific  instructions  are  presented  later.  End  material  data  with 
a  blank  card. 


Projectile  Scale/Shift/Rotate  Card  (7F10.0)  - 


X  SCALE 

Factor  by  which  the  x  coordinates  of  all  projectile  nodes  are  multi¬ 
plied.  Applied  after  the  coordinate  shifts  (XSHIFT,  ZSHIFT) 
described  later.  Must  use  X SCALE  =  YSCALE  =  ZSCALE  =  1.0 
if  any  NABOR  nodes  are  used. 

YSCALE 

= 

Factor  by  which  the  y  coordinates  are  multiplied. 

ZSCALE 

= 

Factor  by  which  the  z  coordinates  are  multiplied. 

XSHIFT 

= 

Increment  added  to  the  x  coordinates  of  all  projectile  nodes 
(length).  Applied  before  the  scale  factors  (XSCALE,  YSCALE, 
ZSCALE). 

ZSHIFT 

= 

Increment  added  to  the  z  coordinates  (length). 

ROTATE 

Rotation  about  the  y  axis  (at  x  =  z  =  0)  of  all  projectile  nodes  (de¬ 
grees).  Applied  after  the  coordinate  shifts  (XSHIFT,  ZSHIFT)  and 
the  scale  factors  (XSCALE,  YSCALE,  ZSCALE).  Clockwise,  as 
viewed  from  a  negative  y  location,  is  a  positive  rotation. 

SLANT 

The  angle  (degrees)  used  to  redefine  the  x  coordinates  of  all  projec¬ 
tile  nodes  with  the  relationship  Xnew  =  X0id  +  Z  tan  (SLANT). 
This  takes  vertical  lines  of  nodes  and  aligns  them  at  an  angle, 
SLANT,  with  the  vertical.  Applied  after  the  other  SC  ALEV 
SHIFT/ROTATE  operations. 

Node  Data  Cards  for  Projectile  -  These  cards  are  required  to  define  the  projectile 
nodes.  If  a  node  is  at  the  interface  of  the  projectile  and  the  target  and  contains  mass  from 
both  the  projectile  and  the  target,  it  must  be  included  with  the  projectile  nodes.  The 
nodes  will  be  numbered  consecutively.  Specific  instructions  for  node  input  data  are 
presented  later.  End  projectile  node  data  with  a  blank  card.  Include  blank  card  even  if 
there  are  no  projectile  nodes. 
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Target  Scale/Shift/Rotate  Card  (7F10.0)  -  Same  as  Projectile  Scale/Shift/Rotate 
Card  except  it  applies  to  the  target  nodes.  Must  be  included  even  if  there  are  no  target 
nodes. 

Node  Data  Cards  for  Target  -  Similar  to  node  data  cards  for  projectile.  Specific  in 
structions  are  presented  later.  End  target  node  data  with  a  blank  card.  Include  blank 
card  even  if  there  are  no  target  nodes. 

Element  Data  Cards  for  Projectile  -  These  cards  are  required  to  define  the  projectile 
elements.  The  element  numbers  will  be  numbered  consecutively.  Specific  instructions 
are  presented  later.  End  projectile  element  data  with  a  blank  card.  Include  blank  card 
even  if  there  are  no  projectile  elements. 

Element  Data  Cards  for  Target  —  Similar  to  element  data  cards  for  projectile. 
Specific  instructions  are  presented  later.  End  target  element  data  with  a  blank  card. 
Include  blank  card  even  if  there  are  no  target  elements. 

MiiHng  Surface  Identification  Card  (1015, 3F10.0)  —  Each  sliding  surface  contains 
one  Identification  Card,  one  or  more  Master  Definition  Cards  and  cards  (if  required) 
describing  the  slave  nodes.  If  there  is  more  than  one  sliding  surface,  all  data  for  the  first 
surface  are  input  before  beginning  data  for  the  second  surface.  The  sliding  surfaces 
should  be  input  in  order  of  the  maximum  node  number  on  each  sliding  surface.  The  mass 
«nrt  sparing  of  the  slave  nodes  should  not  be  significantly  greater  than  that  of  the  master 
nodes  in  the  initial  or  deformed  geometry  unless  a  double  pass  is  used.  Also,  the  slave 
nodes  cannot  be  restrained  in  the  z  direction.  The  nodal  bandwidth  will  be  set  to  contain 
all  sliding  surface  nodes  for  a  specific  sliding  surface.  The  user  should  be  familiar  with 
the  node  generators  before  proceeding. 

MGEOM  =  A  code  describing  the  geometry  of  the  master  surface.  Must  be  1, 

2, 3, 4, 5,  as  shown  in  Figure  17.  Specific  geometries  are  described 
in  the  Master  Definition  Cards. 


'A  code  describing  the  search  routines  used  to  find  the  appropriate 
triangular  plane  on  the  master  surface.  Whenever  possible  use 
the  specialized  searches,  SEEK  =  (+1.  +.2,  ±3).  The  specialized 
routines  can  be  used  whenever  the  master  surface  is  a  single 
valued  function  of  two  coordinates  (i.e.,  any  vertical  line  parallel 
to  a  specified  axis  must  not  pass  through  the  master  surface  at 
more  than  one  point).  For  eroding  sliding  surfaces  set  SEEK  =  4. 

-1  For  all  slave  nodes  on  the  negative  x  side  of  the  master  surface. 

1  For  all  slave  nodes  on  the  positive  x  side  of  the  master  surface. 

-2  For  all  slave  nodes  in  the  negative  y  side  of  the  master  surface. 

=«s  2  For  all  slave  nodes  on  the  positive  y  side  of  the  master  surface. 

-3  For  all  slave  nodes  on  the  negative  z  side  of  the  master  surface. 

3  For  all  slave  nodes  on  the  positive  z  side  of  the  master  surface. 

4  For  the  generalized  search  routine.  Use  of  this  option  can  re¬ 
quire  significant  additional  CPU  time.  Every  master  trian¬ 
gular  surface  is  considered  for  each  slave  node.  If  a  slave  node  is 
contained  in  the  triangular  projection  (onto  a  principal  plane) 
of  one  or  more  master  surface  triangles,  and  if  it  is  close  to  the 
triangular  plane  i  S  n  |  <  REFVEL  -  DT,  then  the  master  plane 
closest  to  the  slave  node  is  selected,  s  n  is  the  normal  distance 
between  the  slave  node  and  the  master  plane,  REFVEL  is 
the  reference  velocity  given  on  this  card  and  DT  is  the  integra¬ 
tion  time  increment.  If  the  slave  node  projection  is  not  within 
any  master  triangular  projections  but  is  close  normally  to  at 
least  one  triangular  plane  |gn  |  <  REFVEL  •  DT,  and  if  the  dis- 
from  the  slave  node  projection  to  the  master  triangular  project- 
tion  is  small,  1 8  edge  I  <  REFVEL •  DT,  then  the  master  triangle 
with  the  smaller  |Sedge|  is  selected.  More  detail  is  given  in 
Reference  10. 

=  Number  of  velocity  iterations.  Errors  in  the  velocity  match  lead  to 
errors  in  the  deviator  and  shear  stresses,  but  generally  not  the 
pressure.  (For  anisotropic  materials,  if  the  pressures  are  deter¬ 
mined  incrementally,  then  there  can  be  errors  introduced  into  the 
pressure.)  For  high  velocity  impact,  where  the  pressures  are 
much  higher  than  the  deviator  and  shear  stresses,  a  relatively 
low  value  of  IT  =  1  or  IT  =  2  can  be  used.  For  lower  pressure 
problems,  higher  values  should  be  used,  IT  =  2  to  IT  =  5.  The 
velocity  iterations,  and  the  corresponding  searches  on  the  master 
surface,  are  performed  only  for  those  slave  nodes  found  to  be  in 
contact  during  the  first  iteration.  For  sliding  surfaces  with  many 
slave  nodes  in  contact  and  many  master  nodes,  high  values  of  IT 


can  lead  to  significant  increases  in  CPU  time.  This  is  especially 
true  if  SEEK  =  4.  For  the  eroding  interface  option  (ERODE  >  0), 
use  IT  =  1. 


Ml 


SI 

SN 

NSN 


NSG 

ISR 


NAB 

FRICT 


=  Node  number  of  the  first  node  on  the  master  surface.  These  will  be 
defined  by  the  Master  Definition  Cards.  The  first  node  is  used  to 
define  the  master  surface,  but  it  may  not  be  the  lowest  numbered 
master  node. 

=  The  lowest  numbered  slave  node.  Leave  blank  for  NSN  =  0. 

=  The  highest  numbered  slave  node.  Leave  blank  for  NSN  =  0. 

=  The  total  number  of  slave  nodes  in  the  sliding  surface.  If  all  nodes 
between  Si  and  SN  are  slave  nodes  (NSN  =  SN  -  Si  + 1),  no  addi¬ 
tional  input  data  are  necessary.  Otherwise,  slave  nodes  are  read 
individually  or  in  groups.  For  NSN  =  0,  the  number  of  slave  nodes 
is  calculated  to  be  the  number  of  nodes  in  the  slave  box  whose 
dimensions  are  given  on  the  Slave  Node  Limits  Cards.  In  some  in¬ 
stances  groups  of  nodes  are  added  to  those  defined  by  the  slave 
box.  Slave  nodes  cannot  be  restrained  in  the  z  direction. 

=  Numbers  of  groups  of  slave  nodes  to  be  read.  If  NSG  =  0,  the  nodes 
are  either  consecutive  between  SI  and  SN,  or  they  are  read  in¬ 
dividually,  or  they  are  designated  by  region. 

=  1  will  release  slave  nodes  from  the  plane  of  symmetry  when  they 
interact  with  master  surface.  May  use  in  conjunction  with 
MGEOM  =  2  for  deep  penetration  problems.  Option  not  used  if 
left  blank. 

=  1  when  the  slave  nodes  are  NAB  OR  nodes;  leave  blank  when 
slave  nodes  are  regular  (not  NABOR)  nodes.  Slave  nodes,  for  a 
specific  sliding  surface,  must  be  all  of  one  type. 

=  The  coefficient  of  sliding  friction. 


E 


REFVEL  =  Reference  velocity,  which  when  multiplied  by  the  integration 

time  increment,  gives  a  reference  distance.  Slave  nodes  are  con¬ 
sidered  to  be  associated  with  a  particular  master  plane  only  when 
they  are  within  this  reference  distance.  It  is  recommended  that 
REFVEL  be  1.0  to  1.5  times  the  relative  impact  velocity,  or  the 
detonation  velocity  of  the  explosive. 

ERODE  =  Equivalent  plastic  strain  (or  volumetric  strain),  which  if  exceeded 

by  any  element  on  the  master  surface,  will  cause  the  element  to  be 
completely  failed.  Subsequently,  the  master  surface  will  be 
redefined  to  go  around  the  failed  element.  This  allows  for  penetra¬ 
tion/perforation  of  thick  plates.  If  ERODE  =  0,  then  erosion  is  not 
used.  The  specific  erosion  algorithm  is  described  in  Reference  8. 
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Because  total  failure  of  the  elements  must  be  achieved  by  the 
eroding  interface  algorithm,  it  is  important  that  EFAIL  (a 
material  property)  be  much  greater  than  ERODE. 

Master  Definition  Card  for  Plate  Geometry  (MGEOM  =  1)  (915)  - 

NML  =  Number  of  nodes  per  row  of  master  nodes.  NML  is  equal  to  N  in 

Figure  17.  Each  row  of  master  nodes  must  have  the  same  number 
of  nodes. 


NMW 


IDL 


IDW 


IDIA 


NBOT 


NLAST 


=  Number  of  rows  of  master  nodes.  NMW  is  equal  to  M  in  Figure  17. 
Note  that  a  properly  described  master  surface  will  pass  the  follow¬ 
ing  test.  Place  a  right-handed  triad  of  orthogonal  vectors  on  node 
Ml.  Point  the  first  vector  (thumb)  away  from  the  master  surface 
towards  the  slave  nodes.  Point  the  second  vector  (index  finger) 
down  the  row  of  nodes  starting  at  Ml.  The  third  vector  (second  fin¬ 
ger)  will  then  point  in  the  direction  of  the  remaining  rows  of 
nodes.  Two  of  the  corner  nodes  will  pass  this  test  when  defined  as 
Ml  and  the  other  two  corner  nodes  will  fail  the  test. 

=  The  node  number  increment  along  the  rows  of  master  nodes.  If 
Ml  =  100,  NML  =  6,  and  IDL  =  2,  then  the  first  row  of  nodes  in 
the  master  surface  consists  of  nodes  100,  102,  104,  106,  108,  110. 
IDL  may  be  negative. 

=  The  node  number  increment  between  the  first  node  of  each  row.  If 
IDW  =  20  and  Ml,  NML  and  IDL  are  as  described  in  the  preceding 
description  of  IDL,  then  the  second  row  of  master  nodes  consists  of 
nodes  120,  122,  124,  126,  128,  130.  IDW  may  be  negative. 

ri  is  for  the  diagonal  orientation  shown  in  Figure  17. 

2  is  for  the  other  diagonal  orientation  where  the  diagonals  go  in 
the  general  direction  from  the  first  master  node  to  the  last  master 
snode. 

=  Lowest  numbered  node  on  the  bottom  surface.  Used  only  for  an 
eroding  plate  to  detect  erosion  through  the  bottom  of  the  plate. 
All  plate  nodes  equal  to  or  larger  than  NBOT  are  considered  to  be 
on  the  bottom  of  the  plate. 

=  The  largest  node  number  which  will  be  on  any  element  on  the  sur¬ 
face  of  an  eroding  sliding  surface.  Required  only  for  an  eroding 
interface.  The  largest  node  number  on  the  plate  is  adequate  but  a 
more  precise  number  can  reduce  the  bandwidth  required  for  the 
sliding  surface.  A  small  estimate  of  NLAST  will  result  in  a  small 
estimate  of  the  required  bandwidth,  but  the  bandwidth  is 
checked  dynamically  during  erosion  and  the  problem  will  stop  if 
not  enough  memory  is  available  for  the  bandwidth  actually  re¬ 
quired. 
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MGNXT  =  MGEOM  for  the  next  geometry  describing  part  of  the  master  sur¬ 

face.  Leave  blank  when  the  master  surface  description  is  com¬ 
plete. 

MlNXT  =  Ml  to  be  used  with  the  next  geometry  specified  by  MGNXT.  Leave 

blank  when  MGNXT  is  blank. 


Master  Definition  Card  for  Nose-Rod  Geometry  (MGEOM  =  2)  (315,  20 X, 
215)  —  This  option  can  be  used  for  deep  penetration  problems  when  the  projectile  is  sig¬ 
nificantly  harder  than  the  target.  The  master  surface  contains  all  external  triangular 
planes  on  the  nose  and  specified  triangles  in  the  rod.  The  nose  and  rod  geometry  must  be 
consistent  with  the  grid  generators  in  Figure  11  and  12.  The  first  master  node,  Ml,  is  at 
the  tip  of  the  nose  as  shown  in  Figure  17. 

For  this  option  the  slave  nodes  are  generally  in  the  target.  It  is  recommended  that  ISR  = 
1  in  the  Sliding  Surface  Identification  Card,  if  there  are  any  restrained  nodes  on  the 
plane  of  symmetry  at  y  =  0.  Use  of  this  option  will  release  the  restraint  when  the  slave 
node  comes  in  contact  with  the  master  surface. 

NOR  =  Outer  node  ring  number  of  the  nose  and  the  rod  of  the  projectile. 

NIR  =  Inner  ring  number. 

NPL  =  Number  of  planes  of  nodes  in  the  rod  included  in  the  master  sur¬ 

face.  The  interface  of  the  rod  and  the  nose  is  considered  to  be 
plane  number  1.  If  NPL  =  2,  then  the  master  surface  would  in¬ 
clude  all  the  triangular  faces  on  the  nose,  plus  those  between  the 
interface  plane  of  nodes  and  the  plane  of  nodes  directly  above  the 
interface  plane. 

MGNXT,  =  Data  for  additional  master  surfaces,  described  previously  for 

MlNXT  MGEOM  =  1. 


Master  Definition  Card  for  Disk  Geometry  (MGEOM  =  3)  (215, 15X,  415)  —  This  op¬ 
tion  can  be  used  if  the  master  surface  is  on  a  disk  whose  nodal  arrangement  is  equivalent 
to  that  of  the  rod  generator.  (A  disk  is  simply  a  very  short  rod  or  cylinder.)  The  first  node, 
Ml,  is  at  the  center  of  the  master  surface  as  shown  in  Figure  17. 

NRING  =  Maximum  node  ring  number  included  in  the  master  surface. 


MCODE 


T  indicates  the  master  surface  is  on  the  top  surface  of  the  disk  (the 
lower  node  numbers  are  on  the  top) 


-1  indicates  the  master  surface  is  on  the  bottom  of  the  disk. 


=  Lowest  numbered  node  on  the  bottom  surface.  Use  only  for  an 
eroding  interface  with  MCODE  =  1  to  detect  erosion  through  the 
bottom  of  the  plate. 


NBOT 


NLAST 


MGNXT, 

M1NXT 


The  largest  node  number  which  will  be  on  the  surface  of  an  erod¬ 
ing  interface.  See  description  given  for  MGEOM  =  1. 

Data  for  additional  master  surfaces,  described  previously  for 
MGEOM  =  1. 


Master  Definition  Card  for  Cylinder  Geometry  (MGEOM  =  4)  (415,  15X,  215)  — 
This  option  can  be  used  if  the  master  surface  is  on  the  inside  of  a  hollow  cylinder  or  the 
outside  of  a  cylinder.  The  nodal  arrangement  of  the  cylinder  must  be  equivalent  to  that 
of  the  rod  generator.  The  first  node,  Ml,  is  on  the  lower  end  of  the  cylinder  as  shown  in 
Figure  17.  (Higher  node  numbers  are  on  the  lower  end.) 

NRING  =  Nodal  ring  number  of  the  surface. 


NPL 

IDL 


MCODE 


=  Number  of  planes  of  nodes  included  in  the  master  surface. 

=  Node  number  increment  in  going  from  one  plane  of  nodes  upward 
to  the  next  plane.  If,  for  instance,  there  are  17  nodes  per  plane, 
then  IDL  =  -17.  IDL  is  always  negative  since  the  increment  goes 
from  the  higher  numbered  nodes  at  the  bottom  to  the  lower  num¬ 
bered  nodes  at  the  top. 


indicates  the  master  surface  is  on  the  inside  of  the  cylinder.  For 
this  option,  the  first  master  node,  Ml,  is  on  the  plane  of  symmetry, 
on  the  positive  x  axis,  when  the  cylinder  is  in  a  vertical  position 
about  the  z  axis. 

-1  indicates  the  master  surface  is  on  the  outside  of  the  cylinder. 
^For  this  option  the  first  node  is  on  the  negative  x  axes. 


MGNXT,  =  Data  for  additional  master  surfaces,  described  previously  for 

M1NXT  MGEOM  =  1. 


Master  Definition  Card  for  General  Geometry  (MGEOM  =  5)  (415)  —  This  option 
can  be  used  when  it  is  necessary  to  describe  a  general  master  surface  which  cannot  be 
defined  by  the  other  master  surface  generators  (MGEOM  =  1,  2, 3, 4).  This  master  sur¬ 
face  is  defined  by  a  series  of  individual  triangular  surfaces,  as  shown  in  Figure  17,  which 
are  defined  in  the  following  card(s).  NBOT,  NLAST,  MGNXT  and  M1NXT  have  all  been 
previously  defined  for  MGEOM  =  1. 

Individual  Master  Surface  Card(s)  for  MGEOM  =  5  (715)  —  This  card  describes  one 
or  more  triangular  surfaces  which  form  all  or  a  part  of  a  master  surface.  Use  as  many 
cards  as  necessary  and  terminate  the  series  of  cards  with  a  blank  card. 
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NCOMP 


=  Number  of  composite  groups  of  triangular  surfaces  to  generate. 
Each  composite  group  contains  one  (M4  =  0),  two  (M4  >  0  and  MS 
=  0),  or  four  (M5  0)  triangles.  In  Figure  17,  NCOMP  =  3  for  both 

cases  shown. 

Ml  =  Number  of  first  node.  Nodes  Ml,  M2,  and  M3  must  be  coun¬ 

terclockwise  when  viewed  from  the  slave  node  side  of  the  master 
surface. 

M2  =  Number  of  second  node. 

M3  =  Number  of  third  node. 

M4  =  Number  of  fourth  node.  If  M4  =  0,  only  one  triangle  will  be 

generated  for  each  composite  group  of  triangles.  If  M4  >  0  and  MS 
=  0,  then  Ml,  M2,  M3  and  M4  must  be  counterclockwise,  and  two 
triangles  (Ml,  M2,  M3  and  Ml,  M3,  M4)  are  generated  for  each 
composite  group,  as  shown  in  Figure  17. 

M5  =  Number  of  fifth  node.  If  M5  =  0,  only  one  or  two  triangles  are 

formed  for  each  composite  group  of  triangles.  If  M5  >  0,  then  four 
triangles  are  generated  for  each  composite  group  of  triangles,  as 
shown  in  Figure  17. 

INC  =  Node  number  increment  (positive  or  negative)  between  cor¬ 

responding  nodes  in  each  composite  group  of  triangular  elements. 
May  be  left  blank  for  consecutive  numbering  (INC  =  1). 

Slave  Node  Limits  Card  for  NSN=0  (6E10.0)  —  All  nodes  found  in  or  on  the  box  are 

taken  to  be  slave  nodes  for  the  sliding  surface.  No  other  slave  node  cards  may  be  used 

with  this  card. 

XMAX  =  Maximum  x  coordinate  of  slave  node  box. 

XMIN  =  Minimum  x  coordinate  of  slave  node  box. 

YMAX  =  Maximum  y  coordinate  of  slave  node  box. 

YMIN  =  Minimum  y  coordinate  of  slave  node  box. 

ZMAX  =  Maximum  z  coordinate  of  slave  node  box. 

ZMIN  =  Minimum  z  coordinate  of  slave  node  box. 


NBOX  =  1  will  read  another  Slave  Node  Limits  Card  to  add  another  box  of 

slave  nodes.  Leave  blank  if  there  are  no  more  boxes  of  slave  nodes 
to  be  added. 

NSG  Grouped  Slave  Node  Cards  for  NSG  >0  (315) 

SlG  =  The  first  node  in  a  group  of  nodes. 

SNG  =  The  last  node  in  a  group  of  nodes. 

INC  =  The  increment  (positive  or  negative)  between  the  slave  nodes.  If 

SlG  =  100,  SNG  =  120,  and  INC  =  5,  then  nodes  100,  105,  110, 
115, 120  will  be  designated  as  slave  nodes. 

Individual  Slave  Node  Cards  for  NSN  >  0  and  NSG  =  0  (1615)  -  Individual  slave 
nodes  are  read  if  the  slave  nodes  are  not  consecutive  between  Si  and  SN,  and  if  they  are 
not  read  in  groups  (NSG  =  0),  and  if  they  are  not  read  in  by  region  (NSN  >  0). 

Si...  SN  =  Individual  slave  nodes  read  in  ascending  order. 

Restrained  Nodes  Identification  Cards  (215, 2X,  31 1)  -  Each  set  of  restrained  nodes 
contains  one  Restrained  Nodes  Identification  Card  and  additional  cards  to  specify  the 
nodes.  The  program  does  not  impose  any  constraint  on  the  number  of  sets  and  each  set 
can  contain  as  many  as  the  node  arrays  can  handle.  If  there  are  no  restrained  node  sets 
(NRST  =  0  in  Miscellaneous  Card),  this  group  of  cards  is  omitted.  If  there  is  more  than 
one  set  of  restrained  nodes,  all  cards  for  the  first  set  are  entered  before  the  Restrained 
Nodes  Identification  Card  for  the  next  set  is  entered.  This  input  redefines  the  restraints 
on  the  designated  nodes  (it  does  not  simply  add  to  existing  restraints).  They  must  be 
read  in  ascending  order. 

NFN  =  Number  of  nodes  in  set. 

NFG  =  Number  of  groups  of  nodes  to  be  read.  If  NFG  =  0  the  nodes  are 

read  individually. 

IX,  IY,  IZ  =1  restrains  nodes  in  X,  Y,  Z  directions,  respectively.  Expanded 

description  given  for  Line  of  Nodes  Description  Card  in  Node 
Geometry  Subsection. 


Individual  Restrained  Nodes  Cards  -  For  NFG  =  0  (1615)  - 

FI. . .  FN  =  Individual  nodes  to  be  restrained  in  ascending  order. 

NFG  Group  Restrained  Nodes  Cards  (315)  - 

FIG  =  First  node  in  the  group  of  nodes  to  be  restrained. 

FNG  =  Last  node  in  the  group  of  nodes  to  be  restrained.  (May  be  zero  if 

there  is  only  one  node  in  the  group.) 

INC  =  Increment  between  nodes  in  the  group  of  restrained  nodes.  (May 

be  zero  if  the  increment  is  one  or  if  no  increment  is  needed.) 

Detonation  Card  (4F10.0)  -  Leave  card  blank  if  no  explosives  are  used. 

XDET  =  x  coordinate  of  detonation  for  high  explosive  (distance). 

YDET  =  y  coordinate  of  detonation  (distance). 

ZDET  =  z  coordinate  of  detonation  (distance). 

TBURN  =  Time  at  which  the  detonation  begins  at  XDET,  YDET,  ZDET. 

Initial  Velocity  Card  (7E10.0)  -  This  card  describes  the  initial  velocity  conditions  of 
the  projectile  and  the  target.  If  there  are  interface  nodes  which  include  mass  from  both 
the  projectile  and  the  target,  the  velocities  of  these  nodes  are  adjusted  by  the  program  to 
conserve  momenta. 

PXDOT  =  Projectile  velocity  in  the  x  direction  (distance/time). 

PYDOT  =  Projectile  velocity  in  the  y  direction. 

PZDOT  =  Projectile  velocity  in  the  z  direction. 

TXDOT  =  Target  velocity  in  the  x  direction. 

TYDOT  =  Target  velocity  in  the  y  direction. 

TZDOT  =  Target  velocity  in  the  z  direction. 

DTI  =  Integration  time  increment  for  the  first  cycle.  This  must  be  less 

than  the  time  required  to  travel  across  the  minimum  altitude  of 
each  tetrahedral  element  at  the  sound  speed  of  the  material  in 
that  element. 


a.  Material  Descriptions 


A  summary  of  input  data  for  material  descriptions  is  given  in  Figure  2.  Four  material 
types  are  available:  solids,  explosives,  crushable  solids  and  anisotropic  solids. 


Material  Cards  for  Solids  from  Library  (415,  2F5.0)  -  Data  for  some  materials  are 
available  from  the  material  library  in  subroutine  MATLIB.  The  specific  materials  are 
listed  as  output  from  the  Preprocessor.  Library  materials  may  be  used  directly  without 
being  called  by  this  card.  If  they  are  called  by  this  card,  however,  the  material  data  will 
be  printed  as  part  of  the  output.  This  card  also  allows  the  user  to  make  some  decisions 
regarding  fracture  options.  The  user  should  read  the  comments  in  subroutine  MATLIB 
to  obtain  the  references  from  which  the  data  were  generated. 


MATL 

0 

IDAM 

IFAIL 

DFRAC 

EFAIL 


=  Material  identification  number.  It  must  be  in  the  range  of  1 
through  50  and  must  correspond  to  a  material  number  in  the 
library. 

=  Code  to  specify  libraiy  material. 


=  1  will  compute  material  damage.  The  option  in  the  library  is 
IDAM  =  0,  which  will  not  compute  damage. 


0  will  not  allow  failure  of  the  material  when  the  damage  exceeds 
1.0,  but  rather  will  continue  to  accumulate  the  damage.  This  is 
the  option  which  is  in  the  library. 

1  will  allow  the  material  to  fracture  partially  when  the  damage 
exceeds  1.0.  Partial  fracture  causes  shear  and  tensile  failure,  so 
vonly  compressive  hydrostatic  pressure  capability  remains. 


=  Factor  by  which  library  fracture  strain  constants  (Dl,  D2, 
EFMIN  -  defined  later)  are  multiplied.  If  left  blank,  (DFRAC  = 
0)  factor  will  be  set  to  DFRAC  =  1.0. 


=  Equivalent  plastic  strain  (or  volumetric  strain),  which,  if  ex¬ 
ceeded,  will  totally  fail  the  element  such  that  it  produces  no 
stresses  or  pressures.  If  left  blank,  the  library  values  of  EFAIL  = 
999  will  govern. 


Five  Material  Cards  for  Solid  Input  Data  [415,  5X,  F5.0,  A48/4(8F10.0)]  —  These  five  cards 
specify  all  the  material  constants  for  a  solid  material.  These  cards  will  supersede  any  material 
library  data  with  the  same  material  number,  MATL.  Only  previously  undefined  variables  will 
be  defined. 


DENSITY 


SPH  HEAT 

CONDUCT 

ALPHA 

TEMPI 

TROOM 

TMELT 

XI 

SHEAR  MOD 

Cl,  C2,  N,  C3, 
M,  C4,  SMAX 


Kl,  K2,  K3 

r 

CL 


=  Material  density  (mass/volume). 

=  Specific  heat  (work/mass/degree). 

=  Thermal  conductivity  (power/distance/degree). 

=  Volumetric  coefficient  of  thermal  expansion  (degree- 1). 

=  Initial  temperature  of  the  material  (degree). 

=  Room  temperature  (degree). 

=  Melting  temperature  of  the  material  (degree). 

=  Extra  material  constant  stored  in  the  C9  array  (not  currently  used). 

=  Shear  modulus  of  elasticity  (force/area). 

=  Constants  to  describe  the  material  strength,  a.The  primary  form  of  the 
strength  equation  is 

a  =  [Cl  +  C2  •  fN][l  +  C3  •  fn<*](l  -  T*M[  +  C4  •  P  (1) 

Where  (  is  the  equivalent  plastic  strain,  f*  =  e/e0  is  the  dimensionless 
strain  rate  for  eQ  =  1.0  s-1,  T*  is  the  homologous  temperature,  and  P  is 
hydrostatic  pressure  (compression  is  positive).  N  must  be  a  positive 
number,  and  the  thermal  softening  fraction,  KT  =  |1  -  T*M],  is  set  to  Kp 
=  1.0  when  M  =  0.  If  SMAX  is  input  as  a  positive  number,  then  the  max¬ 
imum  strength  for  a  is  limited  to  SMAX.  If  left  blank  (SMAX  =  0)  the 
strength  (a)  is  not  limited.  A  more  detailed  description  of  the  strength 
model  and  data  is  given  in  Reference  4. 

A  constant  flow  stress  can  be  obtained  by  setting  Cl  to  the  flow  stress,  N 
=  1.0,  and  C2  =  C3  =  C4  =  SMAX  =  M  =  0.  Cl,  C2  and  SMAX  have 
units  of  stress  (force/area)  and  the  others  are  dimensionless. 

=  Cubic  coefficients  for  the  Mie-Gruneisen  equation  of  state  (force/area) 

P  =  (K1m  +  K2/i2  +  K3/x3X1  -  IV2)  +  Es  (1  +  p)  (2) 

where  p  =  plp0  -  1  and  Es  is  internal  energy  per  initial  volume. 

=  Gruneisen  coefficient  for  Mie-Gruneisen  equation  of  state. 

=  Linear  artificial  viscosity  coefficient  (CL  =  0.2). 


CQ  =  Quadratic  artificial  viscosity  coefficient  (CQ  *  4.0). 

PMIN  =  Maximum  hydrostatic  tension  allowed  (force/area). 

CH  =  Hourglass  artificial  viscosity  coefficient  for  quad  elements  (CH  *  0.02). 

D1...D5  =  Constants  for  the  fracture  model  (Reference  6). 

€f  =  [D1  +  D2expD3<7*|[l  +  D4  •  fm*||l  +  D5  •  T*|  (3) 

Where  t f  is  the  equivalent  strain  to  fracture  under  constant  conditions  of 
the  dimensionless  strain  rate,  «*,  homologous  temperature,  T*,  and  the 
pressure-stress  ratio,  a*  =  am/a.  The  mean  normal  stress  is  am  and  a  is 
the  von  Mises  equivalent  stress.  Expression  is  valid  for  a*  <  1.5.  Dam¬ 
age  is  computed  from  D  =  LAt/e{,  and  fracture  is  allowed  to  occur  when 
D  =  1.0. 

SPALL  =  Tensile  spall  stress  (negative  pressure)  at  which  fracture  can  occur 

(force/area). 

EFMIN  =  Minimum  fracture  strain  allowed.  For  a  *  >  1.5,  varies  linearly  from 

at  a*  =  1.5  to  EFMIN  at  om  =  SPALL 

X2  =  Extra  material  constant  stored  in  the  CIO  array. 

Material  Cards  for  Explosives  from  Library  (215)  —  Similar  to  the  cards  for  the  solid  mate¬ 
rials  in  the  library  except  that  no  options  are  provided  for  fracture. 

Three  Material  Cards  for  Explosives  Input  Data  (215,  20X,  A48/6F10.0/7F10.0)  —  Only 
new  variables  will  be  defined.  See  solid  material  definitions  for  other  variables. 

ENERGY  =  Initial  internal  energy  in  explosive,  E0  (energy/volume). 

DET  VEL  =  Detonation  velocity,  D  (distance/time). 

C1...C5  =  Constants  for  the  JWL  equation  of  state.  If  left  blank,  a  gamma  law 

equation  of  state  is  used  where 

7  =  h  +  D2  p/2E0  (4) 

For  gamma  law  the  pressure  is  determined  from 

P  =  (y-l)E/V  (5) 

For  JWL  the  pressure  is  determined  from 

P  =  Cl  (1-C5/C2  V)  ■  exp(-C2  ■  V) 

+  C3  •  (1-C5/C4  V)  ■  exp(-C4  ■  V) 

+  C5  •  E/V  (6) 

where  E  is  internal  energy  per  initial  volume  and  V  =  V/V0  is  the  rela¬ 
tive  volume.  Cl  ai.d  C3  have  the  units  of  pressure  (force/area)  and  C2, 
C4,  C5  are  dimensionless. 


XI,  X2 


Extra  material  variables  stored  in  arrays  C6  and  C7. 
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Material  Cards  for  Crushable  Solids  from  Library  (216, 15X,  F5.0)  —  These  cards  are  simi¬ 
lar  to  the  cards  for  other  library  materials.  Total  failure  is  allowed  through  EFAIL,  but  frac¬ 
ture  due  to  damage  is  not  allowed. 

Four  Material  Cards  for  Crushable  Solids  Input  Data  [215,  1SX,  F5.0, 
A48/3F10.0/2(8F  1 0.0))  —  These  four  cards  specify  the  material  constants  for  a  crushable  solid 
material  (Reference  5).  Only  new  variables  will  be  defined.  See  previous  material  definitions  for 
other  variables. 

Cl,  C4  =  Constants  to  describe  the  material  strength,  a. 

a  =  Cl  +  C4  •  P  (7) 

SMAX  =  Maximum  strength  allowed  (force/area).  If  left  blank  (SMAX  =  0), 

strength  (a)  is  not  limited. 

Constants  to  describe  the  pressure,  P.  The  model, 
and  specific  data  for  concrete,  are  shown  in  Figure 
18.  The  basic  model  can  also  be  used  for  other 
crushable  solid  materials.  PCRUSH,  Kl,  K2,  K3  and 
KLOCK  have  units  of  pressure  (force/area).  UCRUSH 
and  ULOCK  are  dimensionless. 

XI  =  Extra  material  constant  stored  in  the  Dl  array. 

Eight  Material  Cards  for  Each  Anisotropic  Material  [416,  5X,  F5.0,  A48/3F10.0/2(6F10.0)/ 
2(8F10.0)1  —  Some  special  comments  are  as  follows:  If  Kl  =  0  an  isotropic  pressure,  P*,  is  deter¬ 
mined  from  the  linear  elastic  constants.  For  this  option,  there  are  no  energy  effects  and  the  pres¬ 
sure  is  defined  by 

P  =  P*  +  K2p2  +  K3m3 

where  K2  and  K3  are  isotropic  coefficients  and  p  =  p/p0-l.  Generally,  one  should  set  K2  =  K3  = 
0  since  there  is  no  sound  basis  for  higher  order  isotropic  pressure  terms  in  anisotropic  materi¬ 
als.  If  Kl  =  0  and  r  =  0,  P*  is  computed  directly  from  the  strains  and  is  valid  only  for  small 
strains.  If  Kl  =0  and  T  =  1.0,  P*  is  computed  incrementally  and  large  strains  can  be  accom¬ 
modated  If  incremental  pressures  are  determined  for  elements  containing  nodes  on  a  sliding 
surface,  errors  can  be  introduced  if  the  displacements  are  not  consistent  with  the  velocities. 
This  can  occur  if  there  is  not  a  good  velocity  match  on  the  sliding  surface.  Also,  for  Kl  =  0  and 
T  =  1.0,  no  external  pressures  can  be  applied  since  an  external  pressure  array  is  used  to  store 
P*  between  cycles.  For  Kl  >  0,  the  isotropic  Mie-Gruneisen  Equation  of  State  is  used  for  the 
pressure.  This  option  can  be  used  if  it  is  desirable  to  include  energy  effects  in  the  pressure  com¬ 
putation. 

02.  ^3  =  Euler  angles  between  the  system  axes  and  the  material  axes.  The  materi¬ 

al  input  data  to  follow  should  be  relative  to  the  system  coordinate  axes.  If 
its  initial  orientation  does  not  coincide  with  the  system  axes,  then  the 
Euler  angles  are  used  to  define  the  initial  orientation.  See  Figure  19  for 
definition  of  the  angles.  For  an  anisotropic  rod  which  is  rotated  from  ver¬ 
tical  to  30  degrees  from  vertical,  (0y  -  -30  degrees  in  x-7.  plane)  the  cor 
responding  Euler  angles  are  0 j  =  90°,  0 2  =  -30°  and  0 3  =  -90°. 

vxy,  vyl,  vzx  -  Three  Poisson’s  ratios.  Note  that  vxy  is  defined  as  the  relative  decrease  in 

the  y  direction  due  to  a  tensile  increase  in  the  x  direction. 


PCRUSH, 
UCRUSH, 
Kl,  K2,  K3, 
KLOCK, 
ULOCK 


Ex,  Ey,  E, 

Gxy»  Gy,,  G,x 

X- YIELD 
Y-YIELD 
Z- YIELD 
XY- YIELD 
YZ- YIELD 
ZX- YIELD 


=  Moduli  of  elasticity  in  the  three  principal  directions. 

=  Shear  moduli  of  elasticity  in  the  three  principal  planes  (force/ 
area). 

=  Tensile  flow  stress  in  the  x  direction  (force/area). 

=  Tensile  flow  stress  in  y  direction. 

=  Tensile  flow  stress  in  z  direction. 

=  Shear  flow  stress  in  the  x-y  plane  (force/area). 

=  Shear  flow  stress  in  y-z  plane. 

=  Shear  flow  stress  in  z-x  plane. 


It  should  be  noted  that  the  algorithm  to  find  the  yield  surface  may  not  converge  if  the 
individual  tensile  flow  stresses  (or  individual  shear  stresses)  vary  by  more  than  a  factor 
of  1.5  to  2.5.  When  this  occurs,  an  approximate  radial  return  method  is  used  and  a  warn¬ 
ing  message  is  printed  in  the  cycle  output  data. 

b.  Node  Geometry 


A  summary  of  input  data  for  nodal  geometry  is  given  in  Figure  3.  Nodes  may  be  input  as 
a  line  of  nodes,  and/or  special  shapes  consisting  of  rods,  various  nose  geometries,  flat 
plates  and  spheres.  There  is  no  limit  to  the  number  of  shapes  included  for  the  projectile 
or  the  target.  In  all  cases  the  nodes  are  numbered  consecutively.  If  a  node  is  at  the  inter¬ 
face  of  the  projectile  and  the  target  and  contains  mass  from  both  the  projectile  and  the 
target,  it  must  be  included  with  the  projectile  nodes. 

Line  of  Nodes  Description  Card  (215,  6F10.0,  F7.0,  311)  -  One  card  is  required  for 
each  line  of  nodes  to  be  generated.  Refer  to  Figure  8  for  the  spacing  of  the  nodes. 


1 

= 

Identification  number  for  line  of  nodes  geometry. 

NNODE 

= 

Number  of  nodes  in  the  row  of  nodes.  The  nodes  acre  numbered 
consecutively. 

XI 

= 

x  coordinate  of  first  node  (distance). 

Y1 

= 

y  coordinate  of  first  node  (distance). 

Z1 

at 

z  coordinate  of  first  node  (distance). 

XN 

= 

x  coordinate  of  last  node  (distance).  Leave  blank  if  a  single  node  is 
entered. 

YN 

= 

y  coordinate  of  last  node  (distance). 

ZN 

= 

z  coordinate  of  last  node  (distance). 

EXPAND 

IX 

IY 

IZ 


=  Factor  by  which  the  distance  between  nodes  is  multiplied  going 
from  the  first  node  to  the  last  node.  Leave  blank  for  uniform 
spacing. 

=  1  will  restrain  all  nodes  (N1...NN)  in  the  x  direction.  No  restraint 
if  left  blank. 

=  1  will  restrain  nodes  in  y  direction. 

=  1  will  restrain  nodes  in  z  direction. 


Rod  Node  Identification  Card  (15)  -  The  rod  shape  geometry  descriptions  are  given 
in  Figure  11.  The  rod  is  always  generated  in  a  vertical  position  about  the  z  axis.  When 
viewed  from  the  positive  z  direction,  the  nodes  are  numbered  consecutively  coun¬ 
terclockwise,  inner  to  outer  and  downward.  Only  one-half  the  rod  is  generated  as  shown, 
and  restraints  are  provided  normal  to  the  plane  of  symmetry.  The  rotation  of  the  rod  for 
oblique  impact  is  obtained  with  a  Scale/ShifVRotate  Card. 


2 


Identification  number  for  rod  nodes  geometry. 


Rod  Node  Description  Card  (415,  3F10.0)  - 


NOR 

NIR 

NPLN 

IRAD 

ZTOP 

ZBOT 

EXPAND 


=  Outer  node  ring  number. 

=  Inner  node  ring  number.  The  inner  node  ring  number  for  a  solid 
rod  is  0. 

=  The  number  of  cross-sectional  planes  of  nodes  in  the  rod. 

T)  gives  uniform  radial  spacing  and  constant  z  coordinates  at  the 
top  and  bottom  of  the  rod. 

1  requires  all  r  and  z  coordinates  at  top  and  bottom  of  rod  to  be 
sinput  individually. 

=  The  constant  z  coordinate  of  the  top  of  the  rod  for  IRAD  =  0,  or  the 
top  centerline  z  coordinate  (if  applicable)  for  IRAD  =  1. 

=  The  constant  z  coordinate  at  the  bottom  of  the  rod  for  IRAD  *  0, 
or  the  bottom  centerline  z  coordinate  (if  applicable)  for  IRAD  =  1. 

=  Factor  by  which  the  distance  between  corresponding  nodes  in  the 
vertical  direction  is  multiplied  going  from  top  to  bottom.  Leave 
blank  for  uniform  spacing  in  the  vertical  direction. 


Rod  Node  Radii  Card  for  IRAD  =  0  (4F10.0)  - 
ROTOP  =  Outer  radius  of  the  rod  top. 


RITOP  =  Inner  radius  of  the  rod  top. 
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ROBOT 


Outer  radius  of  the  rod  bottom. 


RIBOT  =  Inner  radius  of  the  rod  bottom. 

Rod  Node  Top  Radii  Card(s)  for  IRAD  =  1  (8F10.0)  - 

RT(I)  =  Top  radii  for  each  node  ring(I  =  NIR,NOR).  Array  RT  is  currently 

dimensioned  for  a  maximum  outer  ring  number  of  16.  IfNIR  =  0, 
then  RT(0)  internally  set  to  0  and  1=1,  NOR. 

Rod  Node  Bottom  Radii  Card(s)  for  IRAD  =  1  (8F10.0)  - 

RB(I)  =  Bottom  radii  for  each  node  ring  (I  =  NIR,  NOR).  Array  RB  is  cur¬ 

rently  dimensioned  for  a  maximum  outer  ring  number  of  16.  If 
NIR  =  0,  then  RB(0)  internally  set  to  0  and  1  =  1,  NOR. 

Rod  Node  Top  Surface  Card(s)  for  IRAD  =  1  (8F10.0)  - 

ZT(I)  =  Top  z  coordinate  for  each  node  ring  (I  =  NIR,  NOR).  IfNIR  =  0 

then  ZT(0)  =  ZTOP  in  the  Rod  Node  Description  Card  and  1  =  1, 
NOR.  Note  that  the  top  z  coordinates  do  not  all  have  to  be  equal  as 
shown  in  Figure  11,  but  can  be  varied  to  form  curved  surfaces. 

Rod  Node  Bottom  Surface  Card(s)  for  IRAD  =  1  (8F10.0)  - 

ZB(I)  =  Bottom  z  coordinate  for  each  node  ring  (I  =  NIR,  NOR).  If  NIR  = 

0,  then  ZB(0)  =  ZBOT  in  the  Rod  Node  Description  Card  and  1  =  1, 
NOR. 

NOTE:  If  it  is  not  possible  to  describe  the  node  geometry  of  the  rod  with  a  single  shape,  it 
is  possible  to  use  multiple  shapes  to  form  a  single  rod.  The  nodes  must  be  numbered 
consecutively,  and  the  radii  and  the  number  of  node  rings  must  be  the  same  for  the  in¬ 
dividual  rod  shapes  at  their  interface.  Also,  ZBOT  and  ZTOP  for  adjoining  rods  should 
not  be  identical.  ZTOP  for  the  lower  rod  should  be  less  than  ZBOT  for  the  upper  rod  by 
the  desired  node  spacing  in  the  z  direction. 

Nose  Node  Identification  Card  (15)  -  The  nose  geometry  descriptions  are  given  in 
Figure  12.  The  nose  shape  is  always  generated  in  a  vertical  position  (pointed  downward) 
about  the  z  axis.  When  viewed  from  the  positive  z  direction,  the  nodes  are  numbered 
consecutively,  counterclockwise,  downward,  and  inner  to  outer.  Only  one-half  the  nose 
shape  is  generated  as  shown,  and  restraints  are  provided  normal  to  the  plane  of  sym¬ 
metry.  The  node  geometry  for  the  plane  of  nodes  at  the  rod  interface  is  not  generated 
with  the  nose  generator  and  must  therefore  be  generated  with  the  rod  generator.  The 
number  of  rings  must  be  identical  for  the  rod  and  the  nose. 

3  =  Identification  number  for  nose  geometry. 

Nose  Node  Description  Card  (415,  4F10.0)  - 

NOR  =  Outer  node  ring  number. 


NIR 


INOSE 


IRAD 

ROTOP 


J 


Inner  node  ring  number. 

1  identifies  a  conical  nose  shape. 


2  identifies  a  rounded  nose  shape.  If  the  length  of  the  nose  is  equal 
to  the  radius,  a  spherical  shape  is  generated.  When  the  length  is 
not  equal  to  the  radius,  the  axial  coordinates  are  scaled  and  the 
various  radii  are  not  changed. 

1^3  identifies  a  tangent  ogival  nose  shape. 


=  Radius  option.  If  IRAD  =  1,  individual  radii  and  z  coordinates  are 
input  for  each  ring  of  nodes.  Leave  blank  (IRAD  =  0)  for  uniform 
spacing. 


=  Top  outer  node  radius  of  nose. 


RITOP  =  Top  inner  node  radius  of  nose. 

ZTOP  =  The  z  coordinate  at  the  top  of  the  nose.  This  is  identical  to  ZBOT 

for  the  rod  interface. 


ZMIN  =  The  z  coordinate  at  the  tip  of  the  nose. 


Nose  Node  Top  Radii  Card(s)  for  IRAD  =  1  (8F10.0)  - 

RT(I)  =  Top  radius  for  each  node  ring  (I  =  NIR,  NOR).  Array  RT  is  cur¬ 

rently  dimensioned  for  a  maximum  outer  ring  number  of  16.  If 
NIR  =  0,  then  RT(0)  internally  set  to  0  and  1  =  1,  NOR. 


Nose  Node  ZMIN  Card(s)  for  IRAD  =  1  (8F10.0)  - 


ZM(I)  =  Minimum  coordinates  for  each  node  ring  (I  =  NIR,  NOR).  Array 

ZM  is  currently  dimensioned  for  a  maximum  outer  ring  number  of 
16.  If  NIR  =  0,  then  ZM(0)  internally  set  to  ZTOP  and  I  =  1,  NOR. 


Flat-Plate  Node  Identification  Card  (215)  -  The  flat-plate  descriptions  are  given  in 
Figure  13.  In  all  cases,  the  lines  connecting  the  a<(jacent  comer  nodes  are  parallel  to  one 
of  the  three  primary  axes.  The  nodes  are  generated  in  rows  parallel  to  the  x  axis  and  are 
numbered  consecutively  within  each  row  in  the  direction  of  the  increasing  x  axis.  The 
TYPE  =  1  option  generates  the  nodes  in  horizontal  planes,  row  by  row  in  the  positive  y 
direction  and  plane  by  plane  in  the  negative  z  direction.  The  TYPE  =  2  option  generates 
the  nodes  in  vertical  planes;  row  by  row  in  the  negative  z  direction  and  plane  by  plane  in 
the  positive  y  direction.  This  TYPE  =  2  option  can  sometimes  be  used  to  minimize  the 
nodal  bandwidth  for  penetration  into  thick  plates. 


4 

TYPE 


=  Identification  number  for  flat-plate  geometry. 
'l  generates  nodes  in  horizontal  planes 
12  generates  nodes  in  vertical  planes 


Flat-Plate  Node  Description  Card  for  TYPE  =  1  (615,  5F10.0)  -  This  option 

generates  nodes  in  horizontal  planes.  The  following  descriptions  refer  to  the  upper  por¬ 
tion  of  Figure  13. 

NX  =  The  number  of  nodes  in  the  x  direction;  in  Figure  13,  NX  =  13. 

NY  =  The  number  of  nodes  in  the  y  direction;  in  Figure  13,  NY  =  7. 

NZ  =  The  number  of  nodes  in  the  z  direction;  in  Figure  13,  NZ  =  5. 

NXEND  =  The  number  of  nodes  in  the  x  direction  in  each  of  the  two  variable 

x  spacing  regions.  The  node  at  the  division  between  the  uniform 
and  the  variable  spacing  sections  is  included  in  this  number.  The 
spacing  is  determined  by  X-EXPAND  and  the  fractional  length 
by  X-PART.  In  Figure  13,  NXEND  =  4.  Depending  on  whether 
NX  is  odd  or  even,  NXEND  can  have  a  maximum  value  of  either 
(NX  +  l)/2  or  NX/2,  respectively,  unless  the  special  option  dis¬ 
cussed  in  X-PART  is  used.  The  remaining  middle  x  region  (if  any) 
is  uniformly  spaced.  Leave  blank  for  uniform  spacing  in  the  x 
direction. 

NYEND  =  The  number  of  nodes  in  the  y  direction  in  the  variable  y  spacing 

region.  The  node  at  the  division  between  the  uniform  and  the 
variable  spacing  sections  is  included  in  this  number.  Spacing  is 
determined  by  Y-EXPAND  and  the  fractional  length  by  Y-PART. 
In  Figure  13,  NYEND  =  4.  NYEND  can  have  a  maximum  value 
of  NY.  The  remainingy  region  (if  any)  is  uniformly  spaced.  Leave 
blank  for  uniform  spacing  in  the  y  direction. 

IY  =1  gives  restraint  in  the  y  direction  only  ify  =  0  (Y1  =  0.0). 

X-EXPAND  =  Factor  by  which  the  x  distance  between  adjacent  nodes  changes 

outward  to  the  ends  for  each  X-PART  variable  spacing  region. 
Same  as  described  for  the  Line  of  Nodes  Description  Card. 

X-PART  =  Fractional  part  of  the  total  x  length  of  the  flat  plate  occupied  by 

each  of  the  variable  x  spacing  regions.  If  X-PART  =  0.0,  the  entire 
spacing  in  the  x  direction  is  uniform.  If  X-PART  =  C.5  and 
NXEND  =  NX/2,  the  spacing  in  the  positive  x  direction  is  vari¬ 
able  for  entire  x  length  (XI  to  XN). 

Y-EXPAND  =  Factor  by  which  the  y  distance  between  adjacent  nodes  changes  in 
the  increasing  y  direction  for  the  Y-PART  variable  spacing 
region. 

Y-PART  =  Fractional  part  of  the  total  y  length  of  the  flat  plate  occupied  by 

the  variable  y  spacing  region. 


Z-EXPAND 


=  Factor  by  which  the  z  distance  between  adjacent  nodes  changes  in 
the  decreasing  z  direction  from  Z1  to  ZN. 


Flat  Plate  Node  Description  Card  for  TYPE  =  2  (615,  5F10.0)  -  This  option 
generates  nodes  in  vertical  planes.  See  lower  portion  of  Figure  13  for  description. 
Definition  of  variables  is  analogous  to  those  of  the  TYPE  =  1  variables. 

Flat-Plate  Node  Size  Card  (6F10.0)  - 


XI 

=  The 

Y1 

=  The 

Z1 

=  The 

XN 

=  The 

YN 

=  The 

ZN 

=  The 

minimum  x  coordinate  of  the  plate  shape, 
minimum  y  coordinate  of  the  plate  shape, 
maximum  z  coordinate  of  the  plate  shape, 
maximum  x  coordinate  of  the  plate  shape, 
maximum  y  coordinate  of  the  plate  shape, 
minimum  z  coordinate  of  the  plate  shape. 


Sphere  Node  Identification  Card  (15)  -  The  cross  section  of  the  bottom  spherical 
shape  is  identical  to  that  shown  for  the  rounded  nose  shape  of  Figure  12.  The  top  one- 
half  cross  section  is  initially  geometrically  symmetric  to  the  bottom.  Only  one-half  of 
the  top  and  bottom  halves  are  generated,  and  restraints  are  provided  normal  to  the  ver¬ 
tical  plane  of  symmetry.  The  sphere  is  generated  with  the  nodes  numbered  as  two 
rounded  circular  noses  having  an  interface  between.  The  top  nose  is  generated  first; 
viewed  from  the  positive  z  direction,  this  generation  is  counterclockwise,  upwards  and 
inner  to  outer.  The  bottom  nose  is  generated  with  the  interface  included  with  each 
spherical  shell;  this  generation  viewed  from  the  positive  z  direction  is  counterclockwise, 
downwards  and  inner  to  outer.  A  summary  of  the  number  of  nodes  included  in  the 
various  shapes  is  given  in  Figure  16. 


5 


Identification  number  for  sphere  geometry. 


Sphere  Node  Description  Card  (215,  5X,  15,  3F10.0)  - 

NOR  =  Outer  node  ring  number. 

NIR  =  Inner  node  ring  number. 


IRAD 

RO 

RI 


"0  gives  RO,  RI,  ZCG  input  option  for  uniform  spacing. 
1  gives  RT(I),  ZCG  input  option  (I  =  NIR,  NOR). 

=  Radius  of  outer  node  ring. 

=  Radius  of  inner  node  ring. 

=  The  z  coordinate  of  the  center  of  the  sphere. 


ZCG 


Sphere  Node  Radii  Card(s)  for  IRAD  =  1  (8F10.0)  - 

RT(I)  =  Radii  for  each  node  ring  (I  =  NIR,  NOR).  Array  RT  is  currently 

dimensioned  for  a  maximum  outer  ring  number  of  16.  If  NIR  =  0, 
then  RT(0)  internally  set  to  0  and  1  =  1,  NOR. 

NABOR  Filled  Plate  Node  Identification  Card  (15)  -  Generates  nodes  for  a  plate 
with  normal  (tetrahedral  element)  edges  and  NABOR  nodes  in  the  center,  as  shown  in 
Figure  14,  so  that  penetration  problems  can  have  NABOR  nodes  in  the  high  distortion 
region  and  normal  elements  in  the  low  distortion  region.  Limited  to  generating  the  +Y 
half  of  mirror  symmetry  about  Y  =  0.  The  node  generation  is  generally  from  the  center 
of  the  plate  to  the  edge.  The  nodes  are  generated  in  a  column  at  a  fixed  X  and  Y  from  the 
largest  Z  to  the  smallest  (or  most  negative)  Z.  The  NABOR  nodes  are  generated  first. 
The  column  at  X  =  (XI  4-  XN)/2,  Y  =  DO/2  is  generated  first  followed  by  the  other  inte¬ 
rior  NABOR  nodes.  The  exterior  layer  of  NABOR  nodes  and  all  the  normal  nodes  are 
generated  in  columns  in  a  generally  counterclockwise  direction.  The  mass  of  the 
NABOR  nodes  is  dependent  only  on  the  density  of  the  material  and  the  diameter  of  the 
NABOR  nodes.  They  cannot  accept  mass  from  attached  elements.  The  diameters  of  all 
NABOR  nodes  in  a  specific  problem  must  be  identical. 

The  size  of  the  NABOR  node  region  can  be  calculated  from  the  number  of  NABOR  nodes 
in  each  direction  and  their  diameter.  See  Figure  14  for  a  typical  example.  The  z  dimen¬ 
sion  has  the  NABOR  nodes  stacked  in  columns  but  alternate  columns  are  offset  by  half  a 
diameter.  The  depth  in  the  z  direction  is  (NNZ+0.5)*DO.  The  y  dimension  consists  of 
layers  of  NABOR  nodes.  The  separation  between  the  layers  is  only  SQRT  (2/3)*DO 
which  makes  the  y  width  ((NNY-l)*SQRT(2/3)+  l)*DO.  The  x  dimension  has  both 
layers  and  an  offset.  The  distance  between  columns  in  the  x  direction  is  SQRT(3/4)*DO. 
The  offset  between  layers  is  (SQRT(l/12))*DO.  The  x  width  is  (NNX*SQRT(3/4)+ 
SQRT(l/12))*DO. 

6  =  Identification  number  for  NABOR  filled  plate  geometry. 

NABOR  Filled  Plate  Node  Description  Card  (715,  5X,  3E10.0)  - 

NNX  =  NumberofNABORnodesinthexdirection.Thisnumbermustbe 

3  or  larger  and  odd.  In  Figure  14,  NNX  =  7. 

NNY  =  Number  of  NABOR  nodes  in  they  direction.  This  number  must  be 

3  or  larger  and  odd.  In  Figure  14,  NNY  =  5. 

NNZ  =  Number  of  NABOR  nodes  in  the  z  direction.  This  number  must  be 

2  or  larger.  In  Figure  14,  NNZ  =  3. 

MATL  =  Material  number  for  plate.  Used  for  NABOR  nodes. 

NXEND  =  Number  of  normal  nodes  on  each  end  in  the  x  direction.  In  Figure 

14,  NXEND  =  3. 

NYEND  =  Number  of  normal  nodes  on  the  positive  y  side  of  the  plate.  Both 

NXEND  and  NYEND  must  be  zero  or  nonzero  together.  In  Figure 
14,  NYEND  =  3. 


IY 


1  gives  restraint  in  the  y  direction  if  Y1  *  0. 


i 

X-EXPAND  =  Factor  by  which  the  *  distance  between  adjacent  nodes  changes 
for  the  regular  (not  NAJBOR)  nodes. 

Y-EXPAND  =  Factor  by  which  the  y  distance  between  adjacent  nodes  changes 
for  the  regular  nodes. 

DO  =  Diameter  of  the  NABOR  nodes.  The  first  DO  read  must  be  posi¬ 

tive.  Following  DOs  may  be  the  same  value  or  the  field  can  be  left 
blank. 

NABOR  Filled  Plate  Size  Card  (5E10.0)  - 
XI  =  The  minimum  x  coordinate  of  the  plate. 

Y1  =  The  minimum  y  coordinate  of  the  plate.  When  NYEND  >  0  then 

Y1  must  be  0.0. 


ZTOP 

XN 


YN 


The  maximum  z  coordinate  of  the  plate. 

The  maximum  x  coordinate  of  the  plate.  The  middle  column  in  the 
plane  of  nodes  which  touches  the  y  =  Y1  plane  is  generated  at 
(X1+XN)*0.5.  When  NXEND  >0,  then  the  width  (XN-Xl)  must 
be  larger  than  (NNX+4)*SQRT(3/4)*DO.  When  NXEND -0,  then 
the  plate  is  all  NABOR  nodes  and  the  width  has  already  been 
determined. 

The  maximum  y  coordinate  of  the  plate.  When  NXEND  >  0,  the 
width  (YN)  must  be  larger  than  (NNY+2)*SQRT(2/3)*DO.  When 
NXEND=0,  then  the  plate  width  has  already  been  determined 
and  YN  is  ignored. 


The  minimum  Z  coordinate  value  is  ZBOT  =  ZTOP  -  (NNZ  +  0.5)*DO. 


NABOR  Brick  Node  Identification  Card  (15)  -  This  option  generates  a  brick  shape 
composed  entirely  of  NABOR  nodes  as  shown  in  Figure  15.  It  generates  the  largest 
NABOR  brick  which  lies  completely  within  the  specified  size.  The  NABOR  nodes  touch 
the  XI,  Yl,  Zl,  planes  but  may  not  reach  to  the  XN,  YN,  ZN  planes. 

7  =  Identification  number  for  NABOR  brick  geometry. 

NABOR  Brick  Node  Description  Card  (215,  E10.0)  - 


MATL  =  Material  number  for  NABOR  nodes. 

IY  =1  gives  restraint  in  the  y  direction  if  Yl  =  0. 


DO  =  Diameter  of  the  NABOR  nodes.  The  first  DO  read  must  be  posi¬ 

tive.  Following  DOs  may  be  the  same  value  or  the  field  can  be  left 
blank. 
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XI 


=  The  minimum  x  coordinate  of  the  brick  shape. 

Y1  =  The  minimum  y  coordinate  of  the  brick  shape. 

Z1  =  The  maximum  z  coordinate  of  the  brick  shape. 

XN  =  The  maximum  x  coordinate  of  the  brick  shape. 

YN  =  The  maximum  y  coordinate  of  the  brick  shape. 

ZN  =  The  minimum  z  coordinate  of  the  brick  shape, 

c.  Element  Geometry 

A  summary  of  input  data  for  element  geometry  is  given  in  Figure  4.  Elements  may  be 
input  as  a  series  of  individual  or  composite  elements  and/or  special  shapes  of  rods,  nose 
geometries,  flat  plates,  and/or  spheres.  The  elements  must  be  assembled  in  a  manner 
consistent  with  the  previously  generated  nodal  geometry.  There  is  no  limit  to  the  num¬ 
ber  of  shapes  included  for  the  projectile  or  the  target. 

Series  of  Composite  Elements  Description  Card  (1215)  -  Element  cards  for  each 
Series  of  Composite  Elements  are  supplied  as  needed  by  the  user.  For  this  discussion  it 
will  be  assumed  that  the  elements  are  entered  as  a  series  of  composite  brick  elements, 
each  containing  six  individual  elements  as  shown  in  Figure  9.  Following  this  im¬ 
mediate  discussion  will  be  an  example  and  instructions  for  generating  a  series  of  in¬ 
dividual  elements. 

1  =  Identification  number  for  a  series  of  composite  elements. 

NCOMP  =  Number  of  composite  elements  in  the  series. 

MATL  =  Material  number  of  the  elements.  If  left  blank,  the  material  num¬ 

ber  from  the  previous  element  data  card  will  be  used. 

N1-N8  =  Node  numbers  of  the  first  composite  brick  element  as  shown  in 

Figure  9.  Nodes  Nl,  N2,  N3,  N4  and  nodes  N5,  N6,  N7,  N8  are 
counterclockwise  when  looking  from  Nl  to  N5. 

INC  *  The  node  number  increment  added  to  the  node  numbers  of  the 

previous  composite  brick  element  for  the  next  composite  brick 
element. 

An  example  of  input  data  for  composite  brick  elements  is  shown  in  Figure  10.  In  the 
upper  left  it  can  be  seen  that  there  are  four  rows  of  nodes  (1  to  4, 5  to  8, 9  to  12, 13  to  16), 
which  are  arranged  to  contain  three  composite  brick  elements.  If  the  first  element  is 
numbered  1,  then  the  first  composite  brick  contains  elements  1  to  6,  the  second  contains 
7  to  12,  and  the  third  contains  13  to  18.  The  first  composite  brick  is  defined  by  nodes  Nl 
=  1 ,  N2  =  5,  N3  =  9,  N4  =  13,  N5  =  2,  N6  =  6,  N7  =  10,  and  N8  =  14.  Note  that  N 1  to  N4 
and  N5  to  N8  are  counterclockwise  when  looking  from  Nl  to  N5.  The  six  individual  ele¬ 
ments  are  generated  according  to  the  arrangement  and  order  (A,  B,  C,  D,  E,  F)  shown  in 


Figure  9.  The  node  numbers  for  each  successive  brick  are  simply  INC  =  1  greater  than 
those  of  previous  brick.  For  the  second  brick,  for  instance,  N1  =  1  +  1  =  2,  N2  =  5+1  = 
6,  N3  =  9  +  1  =  10,  N4  =  13  +  1  =  14,  N5  =  2  +  1  =  3,  N6  =  6+1  =  7,  N7  =  10+1  =  11, 
and  N8  =  14  +  1  =  15. 

It  is  possible  to  generate  a  series  of  individual  tetrahedral  elements  by  letting  N1  to  N4 
be  the  nodes  of  the  first  element,  where  Nl,  N2,  and  N3  are  counterclockwise  when 
viewed  from  N4.  This  option  is  exercised  when  N5  to  N8  are  left  blank.  It  is  also  possible 
to  generate  a  series  of  composite  wedge  elements,  each  containing  three  individual 
tetrahedral  elements.  The  three  elements  in  a  composite  wedge  element  are  numbered 
consecutively.  If  N2  and  N6  are  left  blank,  the  first  three  tetrahedron  elements  (A,  B,  C) 
are  defined  by  nodes  Nl,  N3,  N4,  N5,  N7  and  N8  as  shown  in  Figure  9.  Likewise,  if  N4 
and  N8  are  left  blank,  the  first  three  elements  (D,  E,  F)  are  defined  by  nodes  N 1,  N2,  N3, 
N5,  N6,  and  N7.  Do  not  use  composite  wedge  elements  or  individual  tetrahedral  ele¬ 
ments  in  conjunction  with  the  average  pressure  option  (VRFACT)  described  in  the  Main 
Routine. 

Rod  Element  Identification  Card  (15)  -  Elements  are  generated  for  the  rod  shapes 
illustrated  in  Figure  11  and  described  by  the  Rod  Shape  Node  Cards.  The  elements  are 
numbered  consecutively  and  are  generated  in  layers  of  composite  brick  elements  begin¬ 
ning  with  top  layer  1  and  ending  at  bottom  layer  NLAY.  The  entire  first  layer  of  ele¬ 
ments  is  generated  before  the  second  layer,  etc.,  and  the  composite  brick  elements  of 
each  layer  are  generated  in  a  counterclockwise  manner  for  each  ring  of  elements  from 
the  inner  to  the  outer  ring. 

2  =  Identification  number  for  rod  geometry. 

Rod  Element  Description  Card  (515)  - 
NOER  =  Outer  element  ring  number. 

NIER  =  Inner  element  ring  number.  The  inner  element  ring  number  for  a 

solid  rod  is  1. 

NLAY  =  The  number  of  layers  of  elements  in  the  rod.  The  total  number  of 

elements  in  a  rod  shape  shown  in  Figure  11  is  dependent  on  the 
number  of  layers  and  the  number  of  elements  per  layer.  The  num¬ 
ber  of  elements  per  layer  is  dependent  on  the  inner  and  outer  ele¬ 
ment  ring  numbers.  [For  example:  If  NLAY  =  10,  NIER  =  2  and 
NOER  =  5,  the  total  number  of  elements  can  be  determined 
through  use  of  Figure  16.  The  number  of  elements  per  layer  for 
the  solid  rod  of  NOER  =  5  is  300,  and  of  NIER  =  2  is  48.  Therefore, 
the  total  of  elements  for  the  hollow  cylinder  is  10  (300-48)  =  2520.) 

Nl  =  The  number  of  the  lowest  numbered  rod  node.  For  the  solid  rod, 

this  is  the  centerline  node  on  the  top  end  of  the  rod.  For  the  hollow 
rod,  this  is  the  innermost  clockwise  node  on  the  top  end  of  the  rod 
when  viewed  from  the  top. 

MATL  =  Material  number  of  a  uniform  material  rod.  MATL  =  0  gives  M(I) 

input  option  (I  =  NIER,  NOER),  where  material  numbers  for  each 
element  ring  must  be  input  individually. 


Rod  Element  Material  Card  for  MATL  =  0  (1615)  - 


M(I)  =  Material  number  for  each  element  ring  (I  =  NIER,  NOER).  Array 

M  is  currently  dimensioned  for  a  maximum  outer  ring  number  of 
16. 

Nose  Element  Identification  Card  (15)  —  Elements  are  generated  for  the  nose  shapes 
illustrated  in  Figure  12  and  described  by  the  Nose  Shape  Node  Cards.  The  elements  are 
numbered  consecutively  and  are  generated  in  shells  of  composite  brick  elements  begin¬ 
ning  with  the  innermost  shell  and  ending  with  the  outermost  shell.  The  entire  first  shell 
of  elements  is  generated  before  the  second  shell,  etc.,  and  the  composite  brick  elements 
of  each  shell  are  generated  in  a  counterclockwise  manner  for  each  ring  of  elements  from 
the  top  to  the  bottom  of  each  shell. 

3  =  Identification  number  for  Nose  Geometry. 

Nose  Element  Description  Card  (215,  5X,  215)  - 

NOER  =  Outer  element  ring  number. 

NIER  =  Inner  element  ring  number. 

N 1  =  The  number  of  the  lowest  numbered  nose  node.  (The  nose  does  not 

include  interface  nodes.) 

MATL  =  Material  number  of  a  uniform  material  nose.  MATL  =  0  gives 

M(I)  input  option  (I  =  NIER,  NOER),  where  material  number  for 
each  element  ring  must  be  input  individually. 

Nose  Element  Material  Card  for  MATL  =  0  (1615)  - 

M(I)  =  Material  number  for  each  element  ring  (I  =  NIER,  N OER).  Array 

M  is  currently  dimensioned  for  a  maximum  outer  ring  number  of 
16. 

Flat-Plate  Element  Identification  Card  (215)  -  The  elements  are  generated  for  the 
flat-plate  illustrated  in  Figure  13  and  described  by  the  Flat-Plate  Shape  Node  Cards. 
The  elements  are  numbered  consecutively  and  are  generated  in  rows  of  composite  brick 
elements.  The  rows  of  elements  go  in  the  direction  of  the  increasing  x  axis.  The  TYPE  = 

1  option  generates  the  elements  in  horizontal  layers  and  the  TYPE  =  2  option  generates 
the  elements  in  vertical  layers.  Element  generation  must  be  consistent  with  the  genera¬ 
tion  of  nodes. 

4  =  Identification  number  for  flat-plate  geometry. 


TYPE 


1  generates  elements  in  horizontal  layers. 

2  generates  elements  in  vertical  layers. 


Flat-Plate  Element  Description  Card  (515)  - 

NLX  =  Number  of  layers  of  composite  brick  elements  in  the  x  direction. 

The  total  number  of  nodes  along  the  x  direction  must  be  NLX  +  1. 
In  Figure  13,  NLX  =  12. 


NLY 

NLZ 

N1 


Number  of  layers  of  composite  brick  elements  in  the  y  direction. 
The  total  number  of  nodes  along  the  y  direction  must  be  NLY  +  1. 
In  Figure  13,  NLY  =*  6. 

Number  of  layers  of  composite  brick  elements  in  the  z  direction.  In 
Figure  13,  NLZ  =  4. 

The  node  number  of  the  comer  node  shown  in  Figure  13. 


MATL  =  Material  number  of  the  flat  plate. 


Sphere  Element  Identification  Card  (15)  -  Elements  are  generated  for  a  sphere,  the 
bottom  half  cross  section  of  which  is  identical  to  the  rounded  nose  shown  in  Figure  12 
and  described  by  the  Sphere-Shape  Nose  Cards.  When  viewed  from  the  top,  the  ele¬ 
ments  are  consecutively  numbered  counterclockwise,  upwards  and  outwards  for  the  top 
one-half  and  then  counterclockwise,  downwards  and  outwards  for  the  bottom  one-half. 

5  =  Identification  number  for  sphere  shape. 

Sphere  Element  Description  Card  (215,  5X,  215)  - 

NOER  =  Outer  element  ring  number. 

NIER  =  Inner  element  ring  number.  The  inner  element  ring  number  for  a 

solid  sphere  is  1. 

N1  =  The  number  of  lowest  numbered  sphere  node. 

MATL  =  Material  number  of  a  uniform  material  sphere.  MATL  *  0  gives 

M(I)  input  option  (I  =  NIER,  NOER),  where  material  number  for 
each  element  ring  must  be  input  individually. 

Sphere  Element  Material  Card  for  MATL  =  0(1615)  - 

M(I)  =  Material  number  for  each  element  ring  (I  =  NIER,  NOER).  Array 

M  is  currently  dimensioned  for  a  maximum  outer  ring  number  of 
16. 

NABOR  Filled  Plate  Element  Identification  Card  (15)  -  Elements  are  generated 
for  the  plate  surrounding  the  NABOR  nodes,  as  shown  in  Figure  14.  The  inner  element 
layer  is  generated  between  the  outer  layer  of  NABOR  nodes  and  the  first  inner  layer  of 
normal  nodes.  The  order  of  generation  is  similar  to  the  generation  of  the  normal  nodes, 
in  layers  of  columns  and  generally  counterclockwise. 

=  Identification  for  NABOR  filled  plate  shape. 
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NABOR  Filled  Plate  Element  Description  Card  (715)  - 

NNX  ' 

NNY 

NNZ  >  =  Identical  to  input  used  for  corresponding  NABOR  filled  plate 

MATL  node  description  card. 

NXEND 
NYEND  J 


N1  =  The  number  of  the  lowest  numbered  NABOR  node  in  the  plate. 

2.  INPUT  DATA  FOR  THE  MAIN  ROUTINE 

The  function  of  the  Main  Routine  is  to  perform  the  computations.  The  Main  Routine 
reads  initial  conditions  from  the  restart  tape  which  must  be  generated  from  a  Preproces¬ 
sor  run  or  a  previous  Main  Routine  run.  The  descriptions  that  follow  are  for  the  data  in 
Figure  5.  Consistent  units  must  be  used. 

It  is  possible  to  intellect  user  comments  into  the  data  by  use  of  a  $  character,  as 
described  for  the  Preprocessor. 

Description  Card  (A80)  -  A  description  of  the  problem  provided  by  the  user. 
Identification  Card  (515,  5X,  5F10.0)  - 


CASE 

CYCLE 


IPRES 


CHANG 


IHRG 


=  Case  number  for  run  identification. 

=  The  cycle  number  at  which  the  restart  occurs.  The  cycle  numbers 
for  which  restart  files  are  written  are  given  in  the  printed  output 
of  the  previous  run  (Preprocessor  or  Main  Routine). 

/ 

0  gives  no  applied  pressures  read  or  applied. 

=*s  1  will  use  the  pressure  data  which  was  input  in  a  previous  run. 

2  will  read  applied  pressures  to  be  used  in  subsequent  computa¬ 
tions. 

0  will  not  allow  problem  size  to  be  changed. 

1  will  allow  the  problem  to  be  reduced  in  size  at  a  specified  time. 

=s  2  will  allow  the  problem  size  to  be  reduced  at  one  time,  and  then 
increased  at  another  time.  The  expansions  and  reductions  can  be 
performed  in  either  order.  If  identical  times  are  designated  for  the 
expansion  and  reduction,  the  reduction  will  be  performed  first. 

,3  will  allow  the  problem  size  to  be  expanded  at  a  specific  time. 

=  1  will  use  hourglass  artificial  viscosity  when  computing  brick 
average  pressures  for  solid  materials.  Generally  not  required. 


VFRACT 


«=  Fraction  of  initial  volume  of  a  brick  at  which  individual  element 
pressure  is  computed.  An  average  pressure  is  computed  for  a 
brick  volume  larger  than  VFRACT.  Average  pressures  are  not 
used  if  VFRACT  *0or  the  field  is  blank.  Expanded  description 
given  in  Reference  8.  Use  of  VFRACT  option  requires  all  ele¬ 
ments  to  be  input  with  special  shape  generators  or  composite 
brick  elements.  No  wedges  or  individual  tetrahedral  elements  are 
allowed. 

PMAX  “  Maximum  pressure  allowed  in  any  element.  Pressure  not  limited 

if  PMAX  -  0  (or  blank  field). 

VNREF  *  Reference  closing  velocity  for  NABOR  search  routines.  It  should 

be  0.2  to  1.0  times  the  impact  velocity  (or  relative  closing  velocity) 
in  the  problem.  Lower  values  will  cause  less  searching  but  may 
allow  two  NABOR  nodes  to  become  too  close  before  being 
detected.  Higher  values  will  cause  more  searching  but  will  also 
ensure  that  nodes  do  not  become  too  close  before  being  detected. 

ARMAX  =  For  anisotropic  materials  the  maximum  angle  an  element  can 

rotate  (in  degrees)  before  recalculating  the  element  axes.  A 
properly  chosen  value  can  significantly  reduce  central  CPU  time 
with  only  a  small  loss  of  accuracy.  A  large  value  will  prevent 
recalculating  element  rotations,  which  will  cause  all  elements  to 
be  treated  as  if  they  did  not  rotate.  This  uses  the  minimum  CPU 
time  and  may  be  suitable  for  some  problems.  A  small  value  of  AR¬ 
MAX  -  1.0  degree  will  bypass  the  lengthy  element  axis  calcula¬ 
tion  for  elements  in  regions  where  the  nodes  do  not  have  sig¬ 
nificant  rotations. 

ASERR  *  The  maximum  relative  error  allowed  when  calculating  the  stress 

at  the  yield  surface  for  anisotropic  materials.  The  acceptable  band 
of  stresses  is  from  1.0-ASERR  to  1.0+ ASERR  times  the  actual 
yield  stress.  When  left  blank,  the  value  of  ASERR  =  0.01  will  be 
used. 

Integration  Time  Increment  Card  (6F10.0)  - 

DTMAX  =  The  maximum  integration  time  increment  which  will  be  used  for 

the  equations  of  motion. 

DTMIN  =  The  minimum  integration  time  increment  allowed.  If  exceeded, 

the  results  will  be  written  onto  the  restart  tape,  and  the  run  will 
stop. 

SSF  =  The  fraction  of  the  sound  speed  transit  time  used  for  the  integra¬ 

tion  time  increment.  It  must  be  less  than  unity.  Generally,  SSF  * 
0.9  can  be  used,  but  some  conditions  (such  as  eroding  interfaces) 
may  require  a  lower  value  of  SSF  =  0.5. 


TMAX 


=  The  maximum  time  the  problem  is  allowed  to  run.  This  time 
refers  to  the  dynamic  response  of  the  system,  not  the  central 
processor  time  (CPMAX)  described  next.  The  results  at  time  = 
TMAX  are  written  onto  the  restart  tape,  and  the  run  is  discon¬ 
tinued. 

CPMAX  =  Central  Processor  time  at  which  the  results  will  be  written  onto 

the  restart  tape  and  the  run  will  stop.  This  feature  will  be 
bypassed  if  CPMAX  =  0.  Time  units  are  defined  by  subroutine 
CPCLCK. 

EMAX  =  The  upper  limit  for  total  kinetic  energy  if  applied  pressures  are 

included  (IPRES  =  1  or  2).  This  is  used  for  numerical  instability 
checks.  Run  will  stop  if  the  kinetic  energy  exceeds  EMAX.  Leave 
blank  if  there  are  no  applied  pressures. 

Pressure  Cards  for  IPRES  =  2  (615,  F10.0)  -  These  cards  describe  the  applied  pres¬ 
sures  and  the  elements  to  which  they  are  applied.  If  other  pressures  were  used 
previously,  they  are  all  deleted,  and  the  only  applied  pressures  which  act  are  those  that 
are  input  in  the  current  run.  End  with  a  blank  card. 

ELEl  =  The  first  element  in  a  series  of  elements,  to  which  pressure  is  ap¬ 

plied.  It  must  not  be  less  than  ELEl  or  ELEN  from  a  previous 
Pressure  Card. 

ELEN  =  The  last  element  in  series  of  elements.  It  cannot  be  less  than 

ELEl. 

ELEINC  =  The  element  number  increment  between  ELEl  and  ELEN.  If 

ELEl  =  100,  ELEN  =  120  and  ELEINC  =  5,  then  pressures  are 
applied  to  elements  100,  105,  110,  115,  120. 

N1  =  The  node  number  opposite  the  triangular  face  of  element  ELEl  to 

which  the  pressure  is  applied. 

NN  *  The  node  number  opposite  the  triangular  face  of  element  ELEN 

to  which  the  pressure  is  applied. 

NODINC  =  The  node  number  increment  between  N1  and  NN.  For  the  ele¬ 

ments  described  under  ELE  INC  (100, 105, 110, 115, 120),  if  N1  * 
200,  NN  =  208  and  NODE  INC  =  2,  then  the  pressures  are  ap¬ 
plied  to  the  triangular  faces  opposite  nodes,  200,  202,  204,  206, 
208,  of  elements  100,  105,  110, 115,  120. 

PRES  =  The  pressures  which  are  applied  to  the  triangular  faces  of  the  ele¬ 

ments  described  on  this  card  (force/area). 

Time-Pressure  Cards  for  IPRES  =  2  (2F10.0)  -  These  cards  allow  the  applied  pres¬ 
sures  to  be  varied  as  a  function  of  time.  A  minimum  of  two  cards  must  be  used,  which 


span  the  time  from  the  beginning  of  the  run  to  TMAX.  Program  currently  dimensioned 
for  a  maximum  of  50  cards.  End  with  a  blank  card. 

PTIME  *  The  time  corresponding  to  P(T).  Cards  must  be  input  in  order  of 

increasing  time. 

P(T)  »  The  factor  by  which  all  pressures  are  multiplied  at  the  cor¬ 

responding  time.  Intermediate  values  are  linearly  interpo¬ 
lated  between  values  at  specified  times. 

Change  Card  for  CHANG  =  1  or  2  (F10.0, 1015)  -  This  card  is  used  only  if  changes 
are  to  be  made  which  reduce  the  size  of  the  problem.  The  portions  of  the  problem  which 
remain  are  those  which  were  input  first.  Common  uses  are  to  drop  the  explosive  gases 
after  a  liner  has  been  accelerated,  or  to  drop  the  target  after  a  projectile  has  perforated 
the  target.  For  CHANG  =  1  the  problem  size  is  reduced  only  and  there  are  no  expan¬ 
sions.  For  CHANGE  =  2  there  will  be  problem  expansions  as  well  as  problem  reduc¬ 
tions. 

TDROP  =  Time  at  which  change  occurs. 

NODE  =  Total  number  of  nodes  which  remain  in  the  revised  problem. 

PNODE  =  Number  of  projectile  nodes  which  remain  in  the  revised 

problem. 

ELE  =  Total  number  of  elements  which  remain  in  the  revised 

problem. 

PELE  =  Number  of  projectile  elements  which  remain  in  the  revised 

problem. 

NSLID  =  Number  of  sliding  surfaces  which  remain  in  the  revised 

problem. 

NPLOT  =  Number  of  nodes  for  which  time-history  data  are  written, 

which  remain  in  the  revised  problem. 

LPLOT  =  Number  ofelements  for  which  time-history  data  are  written, 

which  remain  in  the  revised  problem. 

NFAIL  =  Number  of  elements  which  will  be  designated  to  fail  totally. 

This  type  of  failure  sets  all  stresses  in  the  element  to  zero.  It 
essentially  makes  the  element  disappear  except  that  mass  is 
retained  at  the  nodes. 

NAB  =  Total  number  of  NABOR  nodes  which  remain  in  the  revised 

problem. 

PNAB  =  Number  of  projectile  NABOR  nodes  which  remain  in  the 

revised  problem. 


Designated  Element  Failure  Card  for  NFAIL  >  0  (1615)  -  This  card  is  used  only  if 
there  are  elements  to  be  totally  failed. 

EF1 ...  EFN  =  Elements  to  be  totally  failed  in  the  revised  problem. 

Change  Add  Card  for  CHANG  =  2  or  3  (F10.0,6I5)  -  This  card  is  used  only  if  there 
are  additions  to  be  made  to  the  problem.  These  additions  can  be  made  in  conjunction 
with  previous  problem  reductions  (CHANG  =  2)  or  they  can  be  made  to  the  original 
problem  (CHANG  =  3). 

TADD  =  Time  at  which  the  change  (problem  size  expansion)  occurs. 

NSLID  =  Number  of  sliding  surfaces  in  the  expanded  problem.  If  this 

is  greater  than  the  existing  numbers  of  sliding  surfaces,  the 
program  will  read  input  data  for  the  additional  sliding  sur¬ 
faces. 

XYZRIG  =  Code  to  allow  for  rigid  frictionless  surfaces  on  positive  sides 

of  planes  at  x=0,  y=0  and/or  z=0.  See  description  given  in 
the  Miscellaneous  Card  in  the  Preprocessor.  It  must  be  input 
here  even  if  input  in  the  Preprocessor 

SPLIT  =  0  or  1.  See  description  given  in  the  Miscellaneous  Card  in  the 

Preprocessor. 

N1R,  NNR  =  First  and  last  rigid  body  nodes  as  described  in  the  Miscel¬ 

laneous  Card  in  the  Preprocessor.  They  must  be  input  here 
even  if  input  in  the  Preprocessor. 

NIL,  NOR,  NIR,  NPL  =  Node  geometry  descriptions  ifinternal  loads  are  desired.  See 

description  given  in  the  Miscellaneous  Card  in  the 
Preprocessor.  They  must  be  input  here  even  if  input  in  the 
Preprocessor. 

NRST  =  Number  of  groups  of  nodes  to  have  restraints  redefined. 


IPRINT 


0  will  not  print  any  node  or  element  data 
1  will  print  data  for  all  nodes  and  elements 


=<j  2  will  print  data  for  all  nodes  at  y=0  and  all  elements  with 
one  face  aty=0 


3  will  print  data  for  all  added  nodes  and  elements 


4  will  print  data  for  all  added  nodes  aty=0  and  all  added  ele¬ 
ments  with  one  face  at  y=0 


Target  Scale/Shift/Rotate  Card  for  CHANG=2  or  3(F10.0)  -  This  card  is  the  same 
as  used  in  Preprocessor.  Applied  to  added  target  nodes  only  in  the  expanded  problem. 
Must  be  included  even  if  there  are  no  additional  nodes  to  be  input. 
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Additional  Target  Node  Data  for  CHANG =2  or  S  -  These  cards  are  required  to 
define  the  additional  nodes  (for  target  only)  for  the  expanded  problem.  Same  format  as 
used  in  Preprocessor.  End  with  a  blank  card.  Include  blank  card  even  if  there  are  no 
additional  target  nodes. 

Additional  Target  Element  Data  for  CHANG=2  or  3-  These  cards  are  required  to 
define  the  additional  elements  (for  target  only)  for  lthe  expanded  problem.  Same  format 
as  used  in  the  Preprocessor.  End  with  a  blank  card.  Include  blank  card  even  if  there  are 
no  additional  target  elements. 

NSLID,  NRST  data  read  as  required.  Same  format  as  used  in  Preprocessor. 

Target  Velocity  Card  for  CHANG =2  or  3  (SOX,  SF10.0)  -  This  card  defines  the  nodal 
velocities  of  added  target  nodes  only  in  the  expanded  problem.  See  descriptions  for 
Velocity/Detonation  Card  in  Preprocessor.  Must  be  included  even  if  there  are  no  target 
nodes  or  if  no  additional  target  nodes  were  input.  If  there  are  no  target  velocities,  do  not 
leave  the  card  blank.  Instead,  input  velocities  as  0.  Must  be  followed  by  a  blank  card. 

Plot  Card  (515,  5X,  4F10.0)  - 


CYCLE 


LOAD 


NPLOT 


LPLOT 


DT  SYS 


TSYS 


DT  NODE 


Same  restart  cycle  as  specified  on  the  Identification  Card. 
j"0  will  not  write  the  system  data  on  the  Plot  Tape. 

L.1  will  write  all  the  system  data  on  the  Plot  Tape. 

CO  will  not  write  the  internal  loads  data  on  the  Plot  Tape. 

1  will  write  all  the  internal  loads  data  on  the  Plot  Tape. 

This  option  can  only  be  used  if  NIL,  NOR,  NIR  and  NPL  are 
properly  defined  in  the  Preprocessor  Miscellaneous  Card. 

The  number  of  nodes  whose  data  will  be  written  on  the  Plot 
Tape.  The  individual  nodes  are  specified  on  the  Designated 
Nodes  Cards. 

The  number  of  elements  whose  data  will  be  written  on  the 
Plot  Tape.  The  individual  elements  are  specified  on  the 
Designated  Elements  Cards. 

The  time  increment  at  which  the  System  and  Internal  Loads 
data  are  written  on  the  Plot  Tape  (time).  These  quantities  do 
not  vary  as  rapidly  as  do  the  individual  node  and  element 
data  so  a  larger  time  increment  can  be  used. 

The  time  at  which  the  first  System  and  Internal  loads  data 
are  written  on  the  Plot  Tape.  If  left  blank,  the  time  at  the 
beginning  of  the  Main  Routine  run  will  be  used  (time). 

The  time  increment  at  which  the  Individual  Node  and  Ele¬ 
ment  data  are  written  on  the  Plot  Tape  (time).  These  quan- 
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tities  vary  more  rapidly  than  the  System  and  Internal  Loads 
data  so  a  smaller  time  increment  can  be  used. 

TNODE  =  The  time  at  which  the  first  Individual  Node  and  Element 

data  are  written  on  the  Plot  Tape  (time). 

Designated  Nodes  Cards  (1615)  -  This  card  is  used  only  if  there  are  node  data  to  be 
written  on  the  plot  tape  (NPLOT  >  0  on  the  plot  card). 

N1  ...  NN  =  Individual  node  numbers  whose  data  will  be  written  on  the 

Plot  Tape.  They  must  be  in  ascending  order. 

Designated  Elements  Cards  (1615)  -  This  card  is  used  only  if  there  are  element  data 
to  be  written  on  the  plot  tape  (LPLOT  >  0  on  the  plot  card). 

El  ...  EN  =  Individual  element  numbers  whose  data  will  be  written  on 

the  Plot  Tape.  They  must  be  in  ascending  order. 

Data  Output  Cards  (2F10.0,  315)  -  These  cards  are  used  to  specify  various  forms  of 
output  data  at  selected  times,  and  the  last  card  must  be  for  a  time  greater  than  TMAX 
even  though  output  will  not  be  provided  for  that  specific  time.  Recall  that  output  is 
automatically  provided  at  TMAX  and  a  data  output  card  need  not  be  provided  for  this 
time.  End  run  with  a  blank  card. 

TIME  =  Time  at  which  output  will  be  provided. 

ECHECK  =  A  code  which  governs  the  printed  output.  Three  options  are 

provided: 

(1)  If  ECHECK  is  greater  than  1000,  the  individual  node 
and  element  data  will  not  be  printed.  Only  system  data 
such  as  centers  of  gravity,  energies,  momenta,  and  net 
velocities  are  provided  for  the  projectile,  target,  and  total 
system. 

(2)  If  ECHECK  =  999,  system  data  will  be  printed  as  will 
nodal  data  on  the  plane  of  symmetry  at  y=0.  Element 
data  will  not  be  printed. 

(3)  If  ECHECK  is  less  than  999,  the  system  data  and  in¬ 
dividual  node  data  will  be  printed.  Individual  element 
data  will  be  printed  for  all  elements  which  have  an  equiv¬ 
alent  plastic  strain  equal  to  or  greater  than  ECHECK. 
For  example,  if  ECHECK  =  0,  all  element  data  will  be 
printed.  If  ECHECK  =  0.5,  only  those  elements  with 
equivalent  plastic  strains  equal  to  or  greater  than  0.5 
will  have  data  printed. 

NOTE:  There  are  some  instances  when  data  are  printed 
even  though  not  specified  with  a  data  output  card.  If  the  min- 
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LOAD 


SAVE 


YPRINT 


3.  INPUT  DATA  FOR  THE  POSTPROCESSOR 

The  function  of  the  Postprocessor  is  to  provide  plots  of  the  results  in  the  form  of  state 
plots  and  time  plots.  The  state  plots  show  results  for  the  entire  system  at  a  specified  time 
and  the  time  plots  show  results  for  a  specified  variable  as  a  function  of  time. 

It  is  possible  to  interject  user  comments  into  the  data  by  use  of  a  $  character,  as 
described  for  the  Preprocessor. 

a.  State  Plots 

Input  data  for  state  plots  are  summarized  in  Figure  6.  The  first  group  of  data  is  for 
geometry  plots,  velocity  vector  plots  and  contour  plots  of  eight  different  variables. 
These  cards  should  be  input  in  order  of  increasing  cycle  numbers.  This  first  group  of 
cards  (Geometry,  Velocity,  Contour)  is  ended  with  a  blank  card. 

The  second  group  of  cards  is  for  internal  loads  in  a  slender  projectile.  These  should  also 
be  input  in  order  of  increasing  cycle  numbers  and  ended  with  a  blank  card. 

Geometry  Plot  Card  (215,  F10.0,  212,  4X,  212,  18X,  A80)  - 

1  =  Code  to  specify  geometry  plot. 

CYCLE  =  Cycle  number  of  the  plot  which  is  desired.  The  cycle  numbers  of 

the  data  written  on  the  restart  tape  are  given  in  the  printed  out¬ 
put  of  the  Preprocessor  and  the  Main  Routine.  If  CYCLE  -  0,  the 
plots  are  requested  on  the  basis  of  time. 


imum  time  increment,  DTMIN,  is  violated,  or  the  specified 
central  processor  time,  CPMAX,  is  exceeded,  or  the  maxi¬ 
mum  run  time,  TMAX,  is  achieved,  the  results  are  written 
on  the  restart  tape.  The  value  of  ECHECK  used  for  the 
printed  output  is  that  value  which  is  on  the  following  data 
output  card. 

"0  will  not  print  internal  loads  data. 

s  1  will  print  interned  loads  in  a  slender  projectile.  This  option 
can  only  be  used  if  NIL,  NOR,  NIR  and  NPL  are  properly 
defined  in  the  Preprocessor  Miscellaneous  Card. 

f0  will  not  write  results  on  the  restart  tape. 


^1  will  write  results  on  the  restart  tape. 

1  will  restrict  printing  for  node  data  to  nodes  with  y=0,  and 
to  element  data  to  elements  with  one  face  on  the  y =0  plane 
and  an  equivalent  plastic  strain  equal  to  or  greater  than 
ECHECK. 
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TIME 


AXES 


VIEW 


FAIL 

MATL 

TITLE 


=  Time  of  the  plot  which  is  desired.  Plots  can  be  requested  by  either 
TIME  or  CYCLE. 


"0  will  use  the  axes  from  the  previous  plot.  This  option  allows 
deformed  geometry  to  be  plotted  together  with  contours  or 
velocity  vectors,  for  instance.  The  previous  VIEW  is  used. 


1  will  automatically  compute  the  axes  to  include  all  nodes.  The 
vertical  axis  is  specified  to  be  10  units,  and  the  horizontal  axis  is 
13  units,  using  the  same  scale  as  the  vertical  axis.  When 
VIEW  =  5, 6, 7,  then  the  median  of  the  Y,  X,  Z  span  is  used  for  the 
cutting  plane. 

2  will  read  the  coordinate  limits  of  the  plot. 


=  A  code  to  specify  the  view  requested.  VIEW  =  1,  2, 3  requests  two 
dimensional  plots  of  the  x-z,  y-z  and  x-y  axes  respectively.  VIEW 
=  4  requests  a  three-dimensional  plot.  VIEW  *  5,  6,  7  requests 
cutting  plane  plots  with  the  cutting  plane  parallel  to  the  x-z,  y-z, 
and  x-y  planes  respectively.  VIEW  =»  8  requests  a  cutting  plane 
plot  using  an  arbitrarily  positioned  cutting  plane. 

=  1  will  place  an  “x"  at  the  center  of  each  fractured  (partially  failed) 
element  face  for  VIEW  ■  1,  2,  3,  only.  Otherwise  leave  blank. 

will  draw  all  materials. 

(N  will  draw  only  elements  with  material  number  N 
=■  A  title  which  is  printed  on  the  plot. 


Plot  Limits  Card  for  AXES  =  2  (6F10.0)  -  This  card  specifies  the  portion  of  the 
problem  which  is  plotted.  Regions  beyond  those  specified  are  not  plotted.  For  the  two- 
dimensional  plots,  the  vertical  axis  is  10  inches  long,  and  the  horizontal  axis  is  as 
specified.  The  length  of  the  horizontal  axis  will  vary  since  it  has  the  same  scale  as  the 
vertical  axis.  The  vertical  axes  are  the  z  axis  (VIEW  «  1  and  2)  and  the  y  axis  (VIEW  * 
3).  For  the  three-dimensional  plots,  the  axes  are  scaled  such  that  the  entire  region  is 
included  within  the  10-inch  vertical  axis  and  a  13-inch  horizontal  axis. 


XMAX 

XMIN 

YMAX 

YMIN 

ZMAX 


The  maximum  x  coordinate  included  in  the  plot  (length)  When 
VIEW  «  6,  XMAX  is  the  position  of  the  cutting  plane 

The  minimum  x  coordinate  included  in  the  plot 

The  maximum  y  coordinate  included  in  the  plot  When  VIEW  =  5. 
YMAX  is  the  position  of  the  cutting  plane 

The  minimum  y  coordinate  included  in  the  plot 

The  maximum  z  coordinate  included  in  the  plot  When  VIEW  *7. 
ZMAX  is  the  position  of  the  cutting  plane 


ZMIN 

LMAX 


LMIN 


-  The  minimum  i  coordinate  included  in  the  plot. 

-  The  maximum  element  number  included  in  the  plot.  If  both 
LMAX  and  LMIN  are  blank  or  zero,  then  all  elements  are  plotted. 

«  The  minimum  element  number  included  in  the  plot. 


SD  Perspective  Card  for  VIEW  -  4  or  8  (0F1O.0,  IS)  -  This  card  is  included  only  for 
the  3D  plots  and  arbitrarily  oriented  cutting  plane  plots. 

XEYE,  YEYE  *  Coordinates  of  the  observer  (length). 

ZEYE 


XPLANE, 

YPLANE, 

ZPLANE 


Coordinates  included  in  the  plane  on  which  the  results  are  plotted 
for  VIEW  -  4.  The  plane  is  normal  to  a  line  from  XEYE.  YEYE, 
ZEYE  to  XPLANE.  YPLANE.  ZPLANE.  For  VIEW  -  8.  this 
plane  is  the  cutting  plane. 


LHIDE 


[0  will  plot  all  free  surfaces  (no  hidden  lines) 

1  will  plot  only  surfaces  which  are  visible  to  the  observer  This  op 
tion  can  require  significant  CPU  time,  especially  if  the  nodes  are 
not  memory  contained. 

2  will  produce  two  plots,  one  each  with  LHIDE  *  0  and  LHIDE  - 
1 

LHIDE  is  not  used  for  VIEW  -  8 


Velocity  Vector  Plot  Card  (815,  F16.0,  SA9  -  Only  previously  undefined  variables 
are  defined.  See  Geometry  Plot  Card  description  for  others.  Cutting  planes,  VIEW  -  8, 
6,  7,  or  8.  are  not  implemented  for  velocity  plots. 


2 

EDGE 

ARROW 


«  Code  to  specify  velocity  vector  plot. 

✓ 

0  will  draw  all  element  faces. 

-si  will  plot  only  an  outline  around  the  external  surfaces  for 
[VIEW  -  1.  2,  3.  Otherwise  leave  blank. 

»  1  will  place  arrowheads  on  the  velocity  vectors 


VSCALE  -  The  velocity  that  will  give  a  velocity  vector  with  s  length  of  1  0 

using  the  scale  of  the  plot  If  left  blank.  VSCALE  will  be  deter 
mined  to  give  the  longest  vector  a  length  of  two  percent  of  the 
length  of  the  vertical  axis 

Plot  Limit*  Card  for  AXES  =  2  (8F10.0)  -  This  is  the  same  as  described  for  Geometry 
Plots 


SD  Perspective  Card  for  VIEW  «  4  (8P10.0)  -  This  is  *he  same  as  deaenbed  for 

Geometry  Plots  except  hidden  line  option  is  not  allowed 


Contour  Plot  Card  (615, 10X,  5A6, 215)  -  This  card  requests  contour  plots  of  element 
variables.  Contours  are  determined  by  first  computing  the  variable  quantities  at  the 
nodes  (i.e.,  the  nodal  pressure  is  the  average  of  the  pressures  of  all  elements  which  con¬ 
tain  the  node).  Then  the  contours  are  drawn  through  the  nodal  quantities.  The  card  is 
mostly  used  for  two-dimensional  plots,  (VIEW  *  1, 2, 3).  The  3D  plots  have  the  contours 
drawn  only  on  the  element  faces.  Only  previously  undefined  variables  are  defined.  See 
Geometry  and  Velocity  Plot  Cards  for  others.  Cutting  planes,  VIEW  =  5,  6,  7,  or  8,  are 
not  implemented. 


TYPE 

NLINE 


SYMBOL 


PRINT 


A  code  to  specify  which  variable  is  requested.  See  Figure  6  for 
description  of  variables. 

Number  of  contours  to  be  plotted.  Currently  limited  to  8.  If 
NLINE  =  0,  six  contours  will  be  plotted  at  values  of  5,  20,  40,  60, 
80  and  95  percent  of  the  range  between  the  minimum  and  maxi¬ 
mum  variable  quantity  limits. 

Increment  at  which  symbols  are  placed  on  contour  lines. 
SYMBOL  *  I  will  place  symbols  at  the  forward  end  of  each  con¬ 
tour  line  within  an  element,  and  SYMBOL  =  5  will  place  symbols 
at  the  forward  end  of  every  fiflh  element,  etc.  SYMBOL  =  0  will 
not  put  any  symbols  on  the  contour  lines. 

1  will  print  the  nodal  quantities  of  the  specified  variable.  Other¬ 
wise  leave  blank. 


Plot  Limits  Card  for  AXES  =  2  (6F10.0)  -  This  is  the  same  as  specified  for  the 
Geometry  Plot  Card. 

Contour  Specification  Card  (8F10.0)  —  This  card  is  required  only  for  NLINE  >0. 


VAR(I)  -  Magnitude  of  contours  to  be  plotted  (I  =  1,  NLINE). 

Internal  Loads  Plot  Card  (15, 5X,  215, 2X,  SI  1, 15X,  5A6  —  This  card  requests  plots  of 
internal  loads.  The  data  must  have  been  previously  generated  in  a  Main  Routine  run  by 
setting  LOAD  *  I  on  the  Plot  Card.  Sign  conventions  for  the  internal  loads  are  given  in 
Figure  20  The  loads  are  only  for  one-half  the  rod  since  a  plane  of  symmetry  is  used  in 
the  geometry  generators. 
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CYCLE 


Code  which  specifies  internal  loads  plots. 

The  cycle  number  of  the  plot  which  is  desired.  The  cycle  numbers 
of  the  data  written  on  the  plot  tape  are  given  in  the  printed  output 
of  the  Main  Routine.  If  CYCLE  *  0,  the  plots  are  requested  on  the 
basis  of  time. 

Time  of  plot  which  is  desired  Plots  can  be  requested  by  either 
TIME  or  CYCLE 


6.1 


TIME 


AXES 


will  automatically  scale  axes  to  include  all  extreme  data  points. 


i 

! 

I 

I 

I 


IP 

IV 

IM 

TITLE 


^2  will  allow  coordinate  axes  limits  to  be  specified. 

-  1  will  give  plots  of  axial  loads.  Otherwise  leave  blank. 
**  1  will  give  plot  of  shear  loads. 

**  1  will  give  plot  of  bending  moments. 

*  A  title  which  will  be  written  on  the  plot. 


Plot  Limits  Card  for  AXES  =  2  (7F10.0)  - 


PMAX 

PMIN 

VMAX 

VMIN 

BMAX 

BMIM 

LENGTH 


=  Maximum  coordinate  of  vertical  axial  load  axis  (force). 

=  Minimum  coordinate  of  vertical  axial  load  axis. 

«  Maximum  coordinate  of  vertical  shear  load  axis  (force). 

3“  Minimum  coordinate  of  vertical  shear  load  axis. 

*  Maximum  coordinate  of  vertical  bending  moment  axis  (force-dis¬ 
tance). 

=3  Minimum  coordinate  of  vertical  bending  moment  axis. 

-  Maximum  coordinate  of  horizontal  axis,  which  is  defined  as  the 
deformed  centerline  distance  from  the  free  end  of  the  rod  (length). 


b.  Tims  Plots 


Input  data  for  time  plots  are  summarized  in  Figure  7.  System  Plot  Cards  should  be  input 
first,  followed  by  the  Internal  Loads,  Individual  Node  and  Individual  Element  Plot 
Cards.  The  variables  are  plotted  as  a  function  of  time.  The  plot  axes  are  divided  into  ten 
units  each.  End  with  a  blank  card. 


System  Plot  Cards  (215, 5X,  F5.0, 4F10.0,  A20)  -  These  cards  request  plots  of  the  sys¬ 
tem  variables.  Each  plot  contains  data  for  the  projectile,  the  target  and  the  total  system 
(projectile  plus  target).  These  data  must  have  been  previously  written  on  the  plot  tape 
by  setting  SYS  «=  1  on  the  Plot  Card  in  the  main  Routine. 


TYPE 


AXES 


*»  A  code  describing  the  type  of  plot.  See  Figure  7  for  description  of 
type.  It  must  be  in  range  of  1-26. 

(0  will  automatically  select  coordinates  to  include  maximum  and 
minimum  values  of  variable  for  total  duration  of  time. 


SCALE 


1  will  read  the  coordinate  limits  of  the  plot. 

Factor  by  which  the  variables  are  multiplied  before  plotting. 


TMAX 


TMIN 

VMAX 

VMIN 

TITLE 


=  Maximum  time  included  on  horizontal  axis  if  AXES  =  1  (time). 
=  Minimum  time  included  on  horizontal  axis  if  AXES  =  1  (time). 
=  Maximum  variable  included  in  vertical  axis  if  AXES  =  1. 

=  Minimum  variable  included  on  vertical  axis  if  AXES  =  1. 

=  A  title  written  on  the  plot. 


Internal  Loads  Plot  Cards  (315,  F5.0,  4F10.0,  A20)  -  These  cards  request  plots  of 
internal  loads  at  a  specified  location  in  the  slender  projectile.  These  data  must  have 
been  previously  written  on  the  plot  tape  by  setting  LOAD  =  1  on  the  plot  card  in  the 
Main  Routine.  Only  previously  undefined  variables  are  defined.  See  System  Plot  Cards 
description  for  others. 


TYPE 


i 


'27  will  plot  axial  loads. 

28  will  plot  shear  loads. 

29  will  plot  bending  moments. 


LAYER  =  The  layer  oflements  at  which  the  loads  are  plotted.  Layer  I  is  mid¬ 

way  between  node  planes  1-1  and  I  where  node  plane  1  is  at  the 
free  end  of  the  projectile. 

Individual  Node  Plot  Cards  (315,  F5.0,  4F10.0,  A20)  -  These  cards  request  plots  of 
nodal  variables.  These  data  must  have  been  previously  written  on  the  plot  tape  by 
specifying  the  requested  nodes  on  the  Designated  Nodes  Card  in  the  Main  Routine. 
Only  previously  undefined  variables  are  defined. 


TYPE  *  A  code  describing  the  type  of  plot.  See  Figure  7  for  description  of 

types.  It  must  be  in  the  range  of 30-39.  Note  that  acceleration  data 
(TYPE  =  36-38)  may  be  incorrect  for  sliding  surface  and  rigid 
body  nodes. 

NODE  =  Specific  node  for  which  plot  data  are  requested. 

Individual  Element  Plot  Cards  (315,  F5.0, 4F10.0, 3A6)  -  These  cards  request  plots 
of  element  variables.  These  data  must  have  been  previously  written  on  the  plot  tape  by 
specifying  the  requested  elements  on  the  Designated  Elements  Cards  in  the  Main 
Routine.  Only  previously  undefined  variables  are  defined. 

TYPE  *  A  code  describing  the  type  of  plot.  See  Figure  7  for  description  of 

types.  It  must  be  in  the  range  of  40-53. 

ELE  *  The  specified  element  for  which  plot  data  are  requested. 


4.  PROGRAM  STRUCTURE  AND  FILE  DESIGNATION 


EPIC-3  consists  of  four  FORTRAN-77  programs,  PREP,  MAIN,  POST1,  and  P0ST2. 
The  Preprocessor,  PREP,  converts  the  information  given  in  the  input  data  into  the  form 
used  in  a  restart  file  and  writes  a  restart  file.  The  Main  Routine,  MAIN,  reads  a  restart 
file,  advances  the  simulation  for  some  time  cycles,  and  writes  both  a  new  restart  file  and 
plot  files  at  user  requested  simulation  times.  The  Main  Routine  can  read  any  restart  file 
(PREP  or  MAIN),  which  allows  the  simulations  to  be  carried  out  in  several  smaller  runs 
or  one  larger  run.  The  Postprocessor  for  state  plots,  POST1,  produces  plots  which  show 
the  state  of  the  simulation  at  a  given  simulation  time.  POST1  can  read  restart  files  to 
produce  geometry,  velocity  vector,  and  contour  plots.  POST1  can  also  read  the  plot  data 
file  to  produce  internal  loads  plots.  The  Postprocessor  for  time  plots,  POST2,  produces 
plots  which  show  how  simulated  quantities  change  during  the  simulation  time.  POST2 
reads  the  plot  data  file.  The  user  must  request  MAIN  to  write  the  proper  data  into  the 
plot  data  file  so  that  POST2  will  have  the  appropriate  data  available.  The  hierarchy 
charts  for  the  four  programs  are  given  in  Figures  21-24  and  subroutine  groupings  are 
given  in  Tables  1  and  2.  Additional  background  on  program  structure  is  given  in 
References  2  and  9. 


All  the  files  used  by  EPIC-3  are  defined  in  COMMON/FILES/  except  for  two  scratch  files 
used  by  POST1.  The  common  block  is  assigned  values  by  assignment  statements  at  the 
beginning  of  each  program.  The  current  designation  and  descriptions  are: 


IN 

IOUT 

ITAPIN 

ITAPOT 


*  5  Input  file  for  all  programs 
»  6  Output  file  for  all  programs 

=  9  Restart  file  input  read  by  MAIN  and  POST1 

*  10  Restart  file  output  written  by  PREP  and  MAIN 


NFILER 

NFILEW 

LFILER 

LFILER 

ISFILR 

ISFILW 

ITAPLT 

IPLTIN 

ITGAD 


*  1  X  Node  Scratch  file  pair  used  for  out  of  memory 

*  2  >  problems  by  PREP,  MAIN,  POST1 

=  3  X  Element  scratch  file  pair  used  for  out  of 
=  4  ^memory  problems  by  PREP,  MAIN,  POST  1 

=  8  Sliding  surface  scratch  file  used  by  PREP,  MAIN,  POST1 

*  14  Sliding  surface  scratch  file  used  by  MAIN 

*  7  Plot  data  file  written  by  MAIN  and  read  by  POST1  and 

POST2 

■  13  Plot  data  file  read  by  MAIN 

*  11  Target  addition  scratch  file  used  by  MAIN 


Files  1,  2,  3,  4,  8,  11,  and  14  are  opened  as  FORTRAN-77  nameless  scratch  files  when 
needed. 


The  unit  1  is  also  used  by  PREP  and  MAIN.  The  input  file  on  channel  IN  is  stripped  of 
comments  and  written  to  file  TEMP.DAT  on  channel  1.  Then  TEMP.  DAT  is  used  for 
input  on  channel  IN.  TEMP.  DAT  is  deleted  on  normal  program  completion. 

POST1  uses  two  additional  scratch  files. 
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TABLE  1.  SUBROUTINE  GROUPINGS  FOR  EPIC-3  CODE 


PREP 

MAIN1 

MAIN2 

GEN 

NABOR 

POST1 

SUBS 

POST2 

PREP 

MAIN 

MATCH 

BRICK 

CRUSH3 

POST1 

CPCLCK 

POST2 

GEOM 

ASTRES 

MIEGRU 

ECORE 

INVAQ 

ARROW 

ERAFLG 

EPLOT 

MATL 

CHANGE 

MOTION 

EGEOM 

MEGRU3 

AUTOS 

EREAD 

ILPLOT 

MATLIB 

CRUSH 

MOVE 

ENOSE 

NABOR3 

CTOUR 

EREAD1 

NPLOT 

START 

DAMAGE 

MSERA 

ENPLAT 

SEARCH 

CVALUS 

EREAD2 

PDATA 

EGET 

MSTEST 

EPLATE 

STRESS 

DATSET 

ERTAPE 

READTP 

ELOOP 

MTMP2 

EROD 

EDGE 

EWRITE 

SYPLOT 

EPLOTS 

NLOOP 

ESPHER 

GEOM2 

EWTAPE 

EPUT 

NPLOTS 

MASS 

GEOM3 

HDATA 

EQSTS 

ORIENT 

MSFLAG 

GETCON 

LBPOS 

ERODE 

PAVE 

NBCORE 

GPLOT 

LPOS 

EULER 

PCALL 

NBLINC 

HCHECK 

MSNEW 

FIND 

PREAD 

NCORE 

HIDEL 

NBPOS 

FORCE 

REPORT 

NGEOM 

ILOADS 

NFIX 

FSGRAD 

RIGID 

NLINE 

LOADP 

NPOS 

GCON 

ROTOF 

NNBRIK 

PLINE 

NREAD 

HEBURN 

SEEK 

NNOSE 

PRTCON 

NREAD1 

ISTRES 

SLIDE 

NNPLAT 

SETORG 

NREAD2 

LOAD 

STHARD 

NPLATE 

(PERSPC) 

NRTAPE 

STRAIN 

NROO 

SETUP 

NWRITE 

STRESS 

NSPHER 

SORT2 

NWTAPE 

SYMINV 

RGEOM 

SORT3 

RECALL 

SYMOR 

SCYLN 

VPLOT 

SAVE 

TARGAD 

SDISK 

SDATA 

TCHECK 

SETFLG 

TRNST 

SWITCH 

VOLUME 

SLAVEN 

SNNEW 

SNOSE 

SPLATE 

SSEK2 

TETRA 

WEDGE 1 

WEDGE2 

WRrTEG 

N20-010JF 


TABLE  2.  REQUIRED  GROUPS  OF  SUBROUTINES  FOR  VARIOUS  TYPES  OF  RUNS 


TYPE  OF  RUN 


IFILEl 


11 


Used  for  EDGE  option  on  velocity  and  contour  plots 


IFILE2 

=  12; 

ITRAIN 

=  in 

>Used  for  three  dimensional  perspective  plots 

ITEMP 

=  12J 

5.  OPTIONS  FOR  MEMORY  CONTAINMENT  OF  DATA 

EPIC-3  can  simulate  problems  that  are  represented  by  a  number  of  nodes  and  elements 
small  enough  to  fit  into  the  node  and  element  arrays.  But  EPIC-3  can  also  simulate 
problems  whose  representation  exceed  the  capacity  of  the  arrays.  This  capability  is 
achieved  by  using  a  primitive  paging  scheme.  Pages  of  nodes  or  elements  can  be 
shuffled  between  arrays  and  files.  The  price  of  using  this  large  problem  capability  is  a 
slight  increase  in  CPU  time  and  a  large  amount  of  file  access  time.  On  a  typical  modern 
computer  the  file  access  time  for  out  of  memory  elements  roughly  equals  the  CPU  time. 
The  file  access  time  for  out  of  memory  nodes  will  usually  be  less  because  there  usually 
are  fewer  nodes  than  elements.  When  both  nodes  and  elements  are  out  memory,  then 
the  file  accesses  are  synchronized  so  that  the  file  access  times  add  together.  A  notable 
exception  is  in  POST1  when  producing  hidden  line  geometry  plots.  The  hidden  line  algo¬ 
rithm  requires  the  comparison  of  every  node  with  each  exterior  element  face.  A  sig¬ 
nificant  number  of  passes  through  the  node  file  may  be  required  for  one  picture.  The 
most  efficient  way  to  run  a  simulation  is  with  large  enough  array  sites  (nodes  and  ele¬ 
ments)  so  that  the  problem  is  contained  in  memory.  Even  virtual  memory  systems  are 
faster  than  the  EPIC-3  out  of  memory  option.  This  version  of  EPIC-3  has  been  written  to 
allow  array  size  changes  to  be  made  easily. 

The  node  and  element  arrays  in  EPIC-3  are  not  treated  as  single  arrays  but  are  broken 
into  convenient  chunks  called  blocks.  The  paging  scheme  explained  in  the  preceding 
paragraphs  requries  pages  smaller  than  the  array  size.  The  page  size  is  one  block.  The 
vectorization  of  element  calculations  imposes  a  block  size  because  of  the  vectorization 
characteristic  of  computers  and  the  corresponding  temporary  arrays  whch  are  needed. 
The  node  and  element  blocks  can  be  chosen  to  be  different  sizes.  In  both  cases  (node  and 
element)  the  array  sizes  should  be  a  multiple  of  the  corresponding  block  size. 

The  number  of  element  blocks  needed  for  array  containment  is  printed  by  PEEP  after 
the  element  input  data  section.  Since  there  is  no  penalty  for  more  element  blocks  than 
needed  for  array  containment,  a  series  of  problems  can  all  be  run  with  a  version  of 
MAIN  large  enough  to  contain  the  largest  problem.  PREP  executes  quickly  so  there  is 
no  big  advantage  to  array  containment  of  elements  but  it  is  convenient  to  have  PPEP 
arrays  the  same  size  as  MAIN.  PREP  will  never  stop  because  of  problems  with  element 
blocks. 

Choosing  the  number  of  node  blocks  to  use  is  more  complex  when  a  small  node  array 
size  is  desirable.  The  minimum  number  of  node  blocks  required  is  printed  by  PREP  un¬ 
der  the  heading  “Bandwidth  Requirements  for  Main  Routine”.  If  PREP  does  not  con¬ 
tain  this  number  of  node  blocks  it  will  stop  at  that  point.  MAIN  also  requires  the  mini¬ 
mum  number  of  node  blocks. 


Figure  22.  Hierarchy  Chart  for  the  Main  Routine 


ROUTMES 


L 


J 


Figure  24  Hierarchy  Chart  for  the  Time  Plots  Postprocessor 


Eroding  slide  surfaces  pose  a  problem  for  the  smallest  number  of  node  Hocks  required 
If  NLAST,  the  largest  node  number  on  an  eroding  master  surface  at  any  time,  is  known 
then  an  accurate  node  block  number  can  be  given  An  overestimate  ofNLAST  will  over 
estimate  the  node  block  number  If  NLAST  is  not  known  and  zero  given,  then  an  under 
estimate  of  node  blocks  will  be  given  by  PREP,  and  MAIN  will  dynamically  check  as  the 
erosion  happens  When  the  number  of  node  blocks  needed  exceeds  the  number  avail 
able,  then  MAIN  stops.  There  is  no  advantage  to  having  more  than  the  minimum  num 
her  of  node  blocks  until  the  number  sufficient  for  array  containment  is  reached 

6.  INSTRUCTIONS  FOR  CHANGING  PROGRAM  DIMENSIONS 

The  COMMON  declarations  have  been  implemented  in  VAX  sty  le  INCLUDE  files  I  n 
fortunately,  this  is  not  part  of  the  FORTRAN  77  standards,  but  every  mainframe  svs 
tern  has  a  similar  extension  This  method  allows  all  COMMON  declarations  to  be  the 
same  because  they  all  are  copies  of  the  same  file  The  labeled  COMMON  block 
ELEMNT  uses  the  PARAMETER  MAXL  to  declare  the  size  of  the  element  arrays  The 
elements  are  handled  in  blocks  of  size  LBSIZE  Currently  LBSIZE  is  set  to  SO  in  PREP 
Only  complete  blocks  can  be  used  so  a  partial  block  at  the  end  of  arrays  will  l>e  unused 
The  smallest  declaration  for  all  elements  in  memory  is  to  set  MAXL  equal  to  the  smai 
lest  multiple  of  LBSIZE.  which  is  equal  to  or  larger  than  the  number  of  elements  ()ni\ 
MAXL  needs  to  la-  changed  to  change  the  number  of  elements  Alter  MAXL  has  been 
changed,  all  the  programs  except  POST2  should  be  recompiled  ami  linked  Each 
program  will  automatically  switch  to  out  of  memory  node  when  needed  The  element 
array  sizes  in  the  various  programs  iPREP.  MAIN  POST1  d«-  not  need  to  he  four 
d  rated.  but  the  subroutine  packages  C.F.N  and  SUBS  must  hav«  the  same  s;ze  COM 
MON  blocks  as  the  main  program  t  hey  are  twing  1 1  nked  to  The  simplest  t  hi  tig  to  do  is  to 
use  the  same  element  arras  size  everywhere  Similar  considerat  mns  appl\  to  nodes  The 
COMMON  block  is  NODE  with  PARAMETER  MAXN  The  node  block  size  is  NBS1ZE 
which  is  set  to  r>4  m  PREP 


The  block  sizes  have  been  chosen  carefully.  The  CRAY  1,  which  can  vectorize  this  code, 
had  a  limit  of  64  elements  in  a  quick  vector  loop,  so  this  sets  an  upper  limit  on  the  block 
size  The  technique  for  pressure  averaging  in  bricks  requires  that  all  6  elements  be 
present  together  in  a  block,  so  element  blocks  must  be  a  multiple  of  6.  The  existing  ele¬ 
ment  block  size  of  60  is  the  largest  multiple  of  6  equal  to  or  less  than  64.  The  array 
storage  used  to  assist  the  vectorizing  method  in  element  loops  must  be  coordinated  with 
the  element  block  size.  The  PARAMETER  MAXV  used  in  COMMON/VECTOR/  must 
be  equal  to  or  larger  than  LBSIZE.  The  element  and  node  block  sizes,  LBSIZE  and 
NBSIZE.  are  set  in  the  opening  code  of  PREP.  This  is  the  only  place  they  are  set.  They 
are  used  to  define  the  format  of  the  restart  files,  so  they  are  included  in  the  restart  file. 
MAIN  and  POST1  read  LESIZE  and  NBSIZE  from  the  restart  file. 

The  array  storage  used  to  assist  reading  and  writing  node  and  element  blocks  must  also 
be  coordinated  with  block  sizes.  The  element  buffering  arrays  are  contained  in  COM- 
MONEBUF/.  The  minimum  array  sizes  needed  are  PIEB  =  3+(LBSIZE*6)  and 
PREB  =  LBSIZE*25  The  node  buffering  arrays  are  contained  in  COMMON/NBUF/. 
The  minimum  array  sizes  needed  are  PINB  =  2  +  NBSIZE  and  PRNB=NBSIZE*17. 

7  EXAMPLE  PROBLEM 

Preprocessor  and  Main  Routine  input  data  for  a  tungsten  projectile  impacting  a  steel 
target  are  given  in  Figure  25.  The  geometry  of  the  computed  response  is  shown  in 
Figure  26  and  output  data  are  shown  in  Figure  27. 
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Figure  25.  Input  Data  for  the  Example  Problem 
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SECTION  III 

CONCLUSIONS  AND  RECOMMENDATIONS 


User  instructions  have  been  provided  for  the  EPIC-3  code.  This  code  can  be  used  for  a 
wide  range  of  problems  involving  high-velocity  impact  and  explosive-metal  interaction. 
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